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Although the origin of matter-antimatter asymmetry remains unknown, continuing advances in 
theory and improved experimental limits have ruled out some scenarios for baryogenesis, for ex- 
ample the sphaleron baryogenesis at the electroweak phase transition in the standard model. At 
the same time, the success of cosmological inflation and the prospects for discovering supersym- 
metry at the LHC have put some other models in sharper focus. We review the current state of 
our understanding of baryogenesis with the emphasis on those scenarios that we consider most 
plausible. 
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I. INTRODUCTION 

When we observe the universe, the most obvious, and easily studied, objects are stars and gas, made up of protons, 
neutrons and electrons. Astrophysicists speak of the density of protons and neutrons, which constitute the bulk of 
the mass of this matter, as the baryon content of the universe. 
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But we know that there is much more to the universe than baryons. By indirect means, astronomers have estabhshed 
that approximately 1 /3 of the energy density of the universe is in the form of some non-baryonic matter, referred 
to as the dark matter, while 2/3 is in a form with negative pressure, perhaps a cosmological constant. The baryons 
make up a mere 5% of the total energy density of the universe. 

Another, even more striking measure of the smallness of the baryon density is provided by the ratio of baryons to 
photons in the Cosmic Microwave Radiation Background (CMBR). Big Bang nucleosynthesis gives a good measure of 
the baryon density; this measurement is well supported by recent measurements of the fluctuations of the cosmic mi- 
crowave radiation background. As a result, the ratio of baryons to photons is now known to about 5% (,Bennett et al\ . 
120031) : 

^ = (6.1 t°:^)x 10-1". (1) 

There is good evidence that there are no large regions of antimatter at any but cosmic distance 
scales (|Cohen. De Ruiula. and Glashowi Il998^. although some small domains of antimatter in the matter-d ominated 
universe are not ruled out by observations l)Belotskv et al\ . 1 1998': 'Dolg ov and 81113 .119931 iKhloDOv et al ]. l200n>) . 

It was A. Sakharov who first suggested that the baryon density migh t not represent some sort of initial condition, 
but might be understandable in terms of microphysical laws ifSakharovl Il967^ . He listed three ingredients to such an 
understanding: 

1. Baryon number violation must occur in the fundamental laws. At very early times, if baryon number violating 
interactions were in equilibrium, then the universe can be said to have "started" with zero baryon number. 
Starting with zero baryon number, baryon number violating interactions are obviously necessary if the universe 
is to end up with a non-zero asymmetry. As we will see, apart from the philosophical appeal of these ideas, the 
success of inflationary theory suggests that, shortly after the big bang, the baryon number was essentially zero. 

2. CP-violation: If CP (the product of charge conjugation and parity) is conserved, every reaction which produces 
a particle will be accompanied by a reaction which produces its antiparticle at precisely the same rate, so no 
baryon number can be generated. 

3. An Arrow of Time (Departure from Thermal Equilibrium): The universe, for much of its history, was very 
nearly in thermal equilibrium. The spectrum of the CMBR is the most perfect blackbody spectrum measured 
in nature. So the universe was certainly in thermal equilibrium 10^ years after the big bang. The success of 
the theory of big bang nucleosynthesis (BBN) provides strong evidence that the universe was in equilibrium 
two-three minutes after the big bang. But if, through its early history, the universe was in thermal equilibrium, 
then even B and CP violating interactions could not produce a net asymmetry. One way to understand this 
is to recall that the CPT theorem assures strict equality of particle and antiparticle masses, so at thermal 
equilibrium, the densities of particles and antiparticles are equal. More precisely, since B is odd under CPT, its 
thermal average vanishes in an equilibrium situation. This can be generalized by saying that the universe must 
have an arrow of time. 

One of the great successes of the Standard Model is that it explains why baryon and lepton number are conserved, 
to a very good approximation. To understand what this means, consider first the modern understanding of Maxwell's 
equations. A quantum field theory is specified by its field content and by a lagrangian density. In the lagrangian, 
one distinguishes renormalizable and non-renormalizable terms. Renormalizable terms have coefficients with mass 
dimension greater than zero; non-renormalizable terms have coefficients (couplings) with mass dimension less than 
zero. For example, in quantum electrodynamics, the electron mass has dimension one, while the charge of the electron 
is dimensionless (throughout we use conventions where h and c are dimensionless) , so these are renormalizable. In fact, 
requiring Lorentz invariance, gauge invariance, and renormalizability leaves only one possibility for the lagrangian 
of electrodynamics: the Maxwell lagrangian, whose variation yields Maxwell's equations. One can, consistent with 
these symmetry principles, write down an infinite number of possible non-renormalizable terms, which would yield 
non-linear modifications of Maxwell's equations. There is nothing wrong with these, but they are characterized by a 
mass, or inverse length scale, M. So the size of non-linear corrections at wavelength A is of order (AM)~" for some 
integer n. M represents some scale at which the laws of electricity and magnetism might be significantly modified. 
Such corrections actually exist, and are for most purposes quite small. 

Similarly, in the Standard Model, at the level of renormalizable terms, there are simply no interactions one can 
write which violate either baryon number or the conservation of the separate lepton numbers (electron, muon and tau 
number). It is possible to add dimension five operators (suppressed by some scale 1/Af) which violate lepton number, 
and dimension six operators (suppressed by 1/M^), which violate baryon number. Again, these non-renormalizable 
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terms must be associated with some mass scale of some new baryon and lepton violating physics. The dimension five 
lepton-number violating operators would give rise to a mass for the neutrinos. The recent discovery of neutrino mass 
probably amounts to a measurement of some of these lepton number violating operators. The scale of new physics 
associated with these operators can not yet be determined, but theoretical arguments suggest a range of possibilities, 
between about 10" and lO^^ GeV. 

The question, then, is what might be the scale, Mb associated with baryon number violation. At the very least, one 
expects quantum effects in gravity to violate all global quantum numbers (e.g. black holes swallow up any quantum 
numbers not connected with long range fields like the photon and graviton), so Mb < -Mp, where Mp, the Planck 
mass, is about 10^^ GeV, or (2 x 10~'^^cm)~^. The leading operators of this kind, if they have Planck mass coefficients, 
would lead to a proton lifetime of order 10^* years or so. 

If quantum gravitational effects were the only source of baryon number violation, we could imagine that the baryon 
asymmetry of the universe was produced when the temperature of the universe was of order the Planck energy 
(10'^^ °K). Some complex processes associated with very energetic configurations would violate baryon number. These 
need not be in thermal equilibrium (indeed, in a theory of gravity, the notion of equilibrium at such a high temperature 
almost certainly does not make sense) . The expansion of the universe at nearly the moment of the big bang would 
provide an arrow of time. CP is violated already at relatively low energies in the Standard Model (through the 
Kobayashi-Maskawa (KM) mechanism), so there is no reason to believe that it is conserved in very high energy 
processes. So we could answer Sakharov by saying that the magnitude of the baryon number is the result of some 
very complicated, extremely high energy process, to which we will never have experimental access. It might be, in 
effect, an initial condition. 

There are good reasons to believe that this pessimistic picture is not the correct one. First, we are trying to 
understand a small, dimensionless number. But in this Planck scale baryogenesis picture, it is not clear how such a 
small dimensionless number might arise. Second, there is growing evidence that the universe underwent a period of 
infiation early in its history. During this period, the universe expanded rapidly by an enormous factor (at least e^*^). 
Inflation is likely to have taken place well below the scale of quantum gravity, and thus any baryon number produced 
in the Planck era was diluted to a totally negligible level. Third, there are a variety of proposals for new physics - 
as well as some experimental evidence - which suggests that baryon and lepton number violating interactions might 
have been important at scales well below the Planck scale. So there is some reason for optimism that we might be 
able to compute the observed baryon number density from some underlying framework, for which we could provide 
both direct (i.e. astrophysical or cosmological) and/or indirect (discovery of new particles and interactions) evidence. 

Several mechanisms have been proposed to understand the baryon asymmetry: 

1. Planck scale baryogenesis: this is the idea, discussed above, that Planck scale phenomena are responsible for 
the asymmetry. We have already advanced arguments (essentially cosmological) that this is unlikely; we will 
elaborate on them in the next section. 

2. Baryogenesis in Grand Unified Theories (GUT baryogenesis): this, the earliest well-motivated scenario for the 
origin of the asymmetry, will be discussed more thoroughly in the next section. Grand Unified Theories unify the 
gauge interactions of the strong, weak and electromagnetic interactions in a single gauge group. They inevitably 
violate baryon number, and they have heavy particles, with mass of order M^-^^^ w 10^^ GeV, whose decays 
can provide a departure from equilibrium. The main objections to this possibility come from issues associated 
with inflation. While there does not exist a compelling microphysical model for inflation, in most models, the 
temperature of the universe after reheating is well below M^^^. But even if it were very large, there would 
be another problem. Successful unification requires supersymmetry, a hypothetical symmetry between fermions 
and bosons, which will play an important role in this review. Supersymmetry implies that the graviton has a 
spin-3/2 partner, called the gravitino. In most models for supersymmetry breaking, these particles have masses 
of order TeV, and are very long lived. Even though these particles are weakly interacti ng, too nia,ny gr avitinos 
are produced, unless the reheating temperature is well below the unification scale iKallosh et aZ.l 12000(1 . 

3. Electroweak baryogenesis: as we will explain, the Standard Model satisfies all of the conditions for baryogenesis. 
This is somewhat surprising, since at low temperatures the model seems to preserve baryon number, but it turns 
out that baryon and lepton number are badly violated at very high temperatures. The departure from thermal 
equilibrium can arise at the electroweak phase transition ~ a transition between the familiar state in which the 
W and Z bosons are massive and one in which they are massless. This transition can be first order, providing 
an arrow of time. It turns out, however, that as we will explain below, any baryon asymmetry produced is 
far too small to account for observations. In certain extensions of the Standard Model, it is possible to obtain 
an adequate asymmetry, but in most cases the allowed region of parameter space is very small. This is true, 
for example, of the Minimal Supersymmetric Standard Model (MSSM). Experiments will soon either discover 
supersymmetry in this region, or close off this tiny segment of parameter space. 
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4. Leptogenesis: The observation that the weak interactions will convert some lepton number to baryon number 
means that if one produces a large lepton number at some stage, this will be processed into a net baryon and 
lepton number. The observation of neutrino masses makes this idea highly plausible. Many but not all of the 
relevant parameters can be directly measured. 

5. Production by coherent motion of scalar fields (the Affleck-Dine mechanism): This mechanism, which can be 
highly efficient, might well be operative if nature is supersymmetric. In this case, as we will explain in much 
greater detail, the ordinary quarks and leptons are accompanied by scalar quarks and leptons. It has been 
widely conjectured that supersymmetry may be discovered in the next generation of high energy accelerators. 
So again, one might hope to uncover the basic underlying physics, and measure some (but it will turn out not 
all) of the relevant parameters. In non-supersymmetric theories, it is believed that scalar fields with the requisite 
properties (low mass, very flat potentials) are unnatural. This supersymmetric baryogenesis mechanism will be 
the main focus of this review. 

In this review we will survey these, and explain in more detail why the last two are by far the most plausible. 
The question then becomes: can we eventually establish that one or the other is correct? In order to establish 
or rule out particular models for the origin of the matter-antimatter asymmetry, we would hope to bring to bear 
both astrophysical/cosmological and particle physics experiments and observations, as well as theoretical arguments. 
Ideally, we would some day be in the position of measuring all of the parameters relevant to the asymmetry, and 
calculating the asymmetry in much the same way that one presently calculates the light element abundances. One 
question we will ask is: how close can we come to this ideal situation? 

In the next section, after a very brief review of the standard cosmology, we present our survey these mechanisms, 
both explaining how they work and discussing their theoretical plausibility. Both electroweak baryogenesis and 
leptogenesis rely on the existence of processes within the standard model which violate baryon and lepton number at 
high temperatures, and we include a brief explanation of this phenomenon. 

We then turn to a more detailed discussion of coherent production of baryons or leptons, the Affieck-Dine (AD) 
mechanism. This mechanism is potentially extremely efficient; it can also operate relatively late in the history of the 
universe. As a result, it can potentially resolve a number of cosmological puzzles. The AD mechanism presupposes low 
energy supersymmetry. Supersymmetry (sometimes called SUSY for short) is a hypothetical extension of Poincare 
invariance, a symmetry which would relate bosons to fermions. If correct, it predicts that for every boson of the 
standard model, there is a fcrmion, and vice versa. It is believed that the masses of the new particles should be 
about a TeV. As supersymmetry will play an important role in much of our discussion, a brief introduction to 
supersymmetry will be provided in the next section. The supersymmetry hypothesis will be tested over the next 
decade by the Tevatron and the Large Hadron Collider at CERN. Interestingly, most other proposals for baryogenesis 
invoke supersymmetry in some way (including electroweak baryogenesis and most detailed models for leptogenesis). 

II. A BARYOGENESIS ROADMAP 
A. A Cosmology Overview 

Our knowledge of the Big Bang rests on a few key observational elements. First, there is the Hubble expansion 
of the universe. This allows us to follow the evolution of the universe to a few billion years after the Big Bang. 
Second, there is the CMBR. This is a relic of the time, about 10^ years after the Big Bang, when the temperature 
dropped to a fraction of an electron volt and electrons and nuclei joined to form neutral atoms. Third, there is 
the abundance of the light elements. This is a relic of the moment of neutrino decoupling, when the temperature 
was about 1 MeV. As we have noted, theory and observation are now in good agreement, yielding the baryon to 
photon ratio of equation ^ . Finally, there are the fluctuations in the temperature of the microwave background , 
measured recently on angular scales below one degree by BOOMERANG fdc Bcrnardis et a/.', ■2000"Netterfield et all 
l200ai . MAXIMA (S. Hananv et ai, 2000), DASI (Netterfield et ai, 2002), and WMAP (Benn ett et ai . 2003). These 
fluctuations are probably a relic of the era of inflation (discussed in more detail below). The baryon density can 
be inferred independently from the CMBR data and from the BBN determination of the baryon density based on 
the measurements of the primordial deuterium abundance LBurles ,,, NoUett. Turner„ .20 01: Kirkman et ai, 2003). The 
agreement is spectacular: f^sft,^ = 0.0214± 0. 002 based on BBN <|Kirkman et a/.l. l2003j) . while the CMBR anisotropy 
measurements yield flsh^ = 0.0224 ± 0.0009 l|Bennett et alL\2003^ . 

The first and perhaps most striking lesson of the measurements of the CMBR is that the universe, on large scales, 
is extremely homogeneous and isotropic. As a result, it can be described by a Robertson- Walker metric: 




(2) 
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R{t) is known as the scale factor; the Hubble "constant" is H = ^. The Hubble constant is essentially the inverse 
of the time; in the radiation dominated era, H — ^; in the matter dominated era, H — It is puzzling that the 
universe should be homogeneous and isotropic to such a high degree. If one runs the clock backward, one finds that 
vast regions of the universe which have only recently been in cau sal contact have e s sentia l ly the s ame t emperature. 

Inf lation provides an expla i iation for this and other puzz les l|Kolb and Turned . Il99(]l : iLindel . llQQCl) . The basic 
idea recht and SteinhardtL ITosi iGiithl. Il98ll: iLindeL ITosl is that for a brief period, R{t) grew extremely rapidly, 
typically exponentially. This has several effects: 

• The observed universe grew from a microscopically small region, explaining homogeneity and isotropy. 

• k — Q, i.e. the universe is spatially flat. This is now well-verified by observations. 

• Small fluctuations in the metric and the field during inflation explain the observed small (part in 10~^) variation 
in the temperature of the CMBR; detailed features of this structure, in agreement with the inflationary theory, 
have now been observed. These fluctuations provide the seeds for formation of the observed structure in the 
universe. 

• Inflation also explains the absence from the universe of objects such as magnetic monopoles expected in many 
particle physics theories. 

While it is probably fair to say that no compelling microscopic theory of inflation yet exists, as a phenomenological 
theory, inflation is very successful. Most pictures of inflation invoke the dynamics of a scalar fleld in a crucial way. 
This scalar field must have very special properties. Typically, for example, the curvature of its potential must be very 
small. The most plausible theories which achieve this invoke supersymmetry in a significant way. Supersymmetry, 
a hypothetical symmetry between fermions and bosons, will be discussed at greater length later in this article. It 
has been widely considered as a possible solution to many puzzles in particle physics. Most importantly for inflation, 
supersymmetry is a theoretical framework which naturally gives rise to scalars with very flat potentials. It also, almost 
automatically, gives rise to stable particles with just the right properties to constitute the dark matter. There are 
difficulties as well. One is associated with the fermionic partner of the graviton, the gravitino. In many models, this 
particle is very long-lived (r > lO^sec), and can spoil conventional Big Bang cosmology if too many are produced. 

There is not space here to review the subject of inflation. Instead, we will give a "narrative" of a possible history 
of the universe, which will be useful to orient our discussion: 

• Before t « lO^'^^ seconds, the universe was very inhomogeneous, with an extremely large energy density. At 
t K, 10~^^, inflation began in a small patch. This was associated with a scalar field, called the inflaton, which 
moved slowly toward the minimum of its potential. 

• The scale of the infiaton potential was of order I0^° GeV^, give or take a few orders of magnitude. 

• During infiation, the scale factor increased by an enormous factor. Any conserved or approximately conserved 
charges, such as monopole number or baryon number, were reduced by at least a factor of 10^'' in this process 

• Infiation ended as the inflaton approached the minimum of its potential. At this point, decays of the inflaton 
lead to reheating of the universe to a high temperature. Depending on the detailed microscopic picture, there are 
constraints on the reheating temperature. If nature is supersymmetric, there is often a danger of producing too 
many gravitinos and other long-lived particles. Typically, this constrains the reheating temperature to be below 
10^ GeV. Even without supersymmetry, detailed inflationary models have difficulty producing high reheating 
temperatures without fine tuning. 

• The baryon asymmetry is generated some time after the era of infiation. Any upper limit on the reheating 
temperature constrains the possible mechanisms for baryogenesis. 



B. Planck Scale Baryogenesis 

It is generally believed that a quantum theory of gravity cannot preserve any global quantum numbers. For example, 
in the collapse of a star to form a black hole, the baryon number of the star is lost; black holes are completely 
characterized by their mass, charge and angular momentum. Virtual processes involving black holes, then, would also 
be expected to violate baryon number. 

In string theory, the only consistent quantum the ory of gravity we k now, these prejudices are born out. There 
are no conserved global symmetries in string theorv ljBanks et o/l llQSSfl . While we can't reliably extract detailed 
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predictions from quantum gravity for baryon number violation, we might expect that it will be described at low 
energies by operators which appear in an effective field theory. The leading operators permitted by the symmetries 
of the Standard model which violate baryon number carry dimension six. An example is: 

Cb = ■^^d*d*d* (3) 

In this equation and those which follow, the various fermion fields, d d, e, e, etc. are spinors of left-handed chirality. 
d contains the creation operator for the right-handed d quark; d^ for the left-handed anti-d quark. The other two 
d-quark states are created by d and . Because the operator is of dimension six, we have indicated that its coefficient 
has dimensions of inverse mass-squared. This is analogous to the effective interaction in the Fermi theory of weak 
interactions. If quantum gravity is responsible for this term, we might expect its coefficient to be of order l/M^, 

where Mp = y^G^ = IQi^GeV. 

Because of this very tiny coefficient, these effects could be important only at extremely early times in the universe, 
when, for example, H ^ Mp. It is probably very difficult to analyze baryon production in this era. It is certainly 
unclear in such a picture where the small number 10"^" might come from. But even if the baryon number was produced 
in this era, it was completely washed out in the subsequent period of inflation. So gravitational baryogenesis seems 
unlikely to be the source of the observed matter-antimatter asymmetry. 



C. GUT baryogenesis 

The earliest well-motivated scenarios for i mplementing Sakharov' s ideas within a detailed microscopic theory were 
provided by grand unified theories ('GUTs') ijKolb and Turneil Il99(]() . In the Standard Model, the strong, weak and 
electromagnetic interactions are described by non-abefian gauge theories based on the groups SU{3), SU{2) and C/(l). 
Grand unification posits that the underlying theory is a gauge theory with a simple group; this gauge symmetry is 
broken at some very high energy scale down to the group of the Standard Model. This hypothesis immediately 
provides an explanation of the quantization of electric charge. It predicts that, at very high energies, the strong, weak 
and electromagnetic couplings (suitably normalized) should have equal strength. And most important, from the point 
of view of this article, it predicts violation of baryon and Icpton numbers. 

If nature is not supcrsymmctric, the GUT hypothesis fails. One can use the renormalization group to determine 
the values of the three gauge couplings as a function of energy, starting with their measured values. One finds that 
they do not meet at a point, i.e. there is no scale where the couplings are equal. Alternatively, one can take the best 
measured couplings, the SU{2) and f/(l) couplings, and use the G UT hypothesis to predict th e value of the strong 
coupling. The resulting prediction is off by 12 standard deviations ijParticle Data Groud.l20n2t) . But if one assumes 
that nature is supersymmetric, and that the new particles predicted by supersymmetry all have masses equal to 1 TeV, 
one obtains unification, within 3 a. The scale of unification turns out to be M^-^^,^ « 2 x 10"'^^. Relaxing the assumption 
that the new particles are degenerate, or assuming that there are additional, so-called threshold corrections to the 
couplings at the GUT scale (~ 4%) can yield complete agreement. 

This value of M^^^ is quite interesting. It is sufficiently below the Planck scale that one might hope to analyze these 
theories without worrying about quantum gravity corrections. Moreover, it leads to proton decay at a rate which may 
be accessible to current proton decay experiments. In fact, the simplest SUSY GUT, based on the ga uge group SU{h) 
is almost completely ruled out by the recent Super-Kamiokande bounds jMuravama and Piercell2002|) . However, there 
are many other models. For example, non-mini mal SU(5) or SO(IO) SUSY GUTs may have a proton lifetime about a 
factor of 5 above the presen t experimental limit j Altarelli^ Feruglio and Masinall200dlBabu. Pati and Wilczekll200d 
iBaic. Perez and Se nianovi?. '2002*: 'Dermisek.Mafi and Rab yl l200lf) . Witten has recently advocated an approach to 
GUT model building (Fricdmann and Witten, 2002; Wittenl . [2002|) which resolves certain problems with these models, 
and in which proton decay might be difficult to see even in large detectors which are being considered for the future. 

GUTs provide a framework which satisfies all three of Sakharov's conditions. Baryon number violation is a hallmark 
of these theories: they typically contain gauge bosons and other fields which mediate B violating interactions such as 
proton decay. CP violation is inevitable; necessarily, any model contains at least the KM mechanism for violating CP, 
and typically there are many new couplings which can violate CP. Departure from equilibrium is associated with the 
dynamics of the massive, B- violating fields. Typically one assumes that these fields are in equilibrium at temperatures 
well above the grand unification scale. As the temperature becomes comparable to their mass, the production rates 
of these particles fall below their rates of decay. Careful calculations in these models often leads to baryon densities 
compatible with what we observe. 

We can illustrate the basic ideas with the simplest GUT model, due to iGeorgi and Glashowl l)l974() . Here the 
unifying gauge group is SU{b). The model we will discuss is not supersymmetric, but it illustrates the important 
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features of GUT baryon number production. A single generation of the standard model (e.g. electron, electron 
neutrino, u-quark, and d-quark) can be embedded in the 5 and 10 representation of SU{5). It is conventional, and 
convenient, to treat all quarks and leptons as left-handed fields. So in a single generation of quarks and Icptons one 
has the quark doublet, Q, the singlet u and d antiquarks {their antiparticles are the right-handed quarks), and the 
e 



lepton doublet, L = 



Then it is natural to identify the fields in the 5 as 



d ' 

d 
e 



(4) 



The generators of SU{3) of color are identified as: 



T : 



(5) 



where A" are the Gell-Mann matrices, while those of SU{2) are identified with: 



T = 




4 



(6) 



The U{1) generator is 



Y' 



/2 



V 



-3 



(7) 



-3/ 



Here the coefficient has been chosen so that the normalization is the same as that of the SU (3) and SU{2) matrices 
(Tr(T°T'') = Sab)- The corresponding gauge boson couples with the same coupling constant as the gluons and W and 
Z bosons. This statement holds at M^^^ ; at lower energies, there are significant radiative corrections (which in the 
supersymmetric case reproduce the observed low energy gauge couplings). 

In the standard model, the hypercharge is related to the ordinary electric charge, Q, and the isospin generator, 
T3, by Q = T3 + y. So one sees that electric charge is quantized, and that Y = ■\/3/40 Y' . Since Y couples 
with the same strength as the SU{2) generators, this gives a prediction of the U{1) coupling of the standard model, 
and correspondingly of the Weinberg angle, sm'^{6w) = 3/8. This prediction receives radiative corrections, which, 
assuming supersymmetry, bring it within experimental errors of the measured value. 

In a single generation, the remaining fields lie in the 10 representation. The 10 transforms as the antisymmetric 
product of two 5's. It has the form 



10,; 



f 





U2 


-Ul 


Q\ 


Ql\ 




-U2 







Ql 


Ql 




Ul 







Ql 


Ql 




-Q\ 


-Qh 


-Ql 





e 


{ 


-Ql 


-Ql 


-Ql 


— e 


/ 



(8) 



where = (?i*, d') are left-handed quark fields, which transform as doublets under the SU(2). SU{5) is not a manifest 
symmetry of nature. It can be broken by the expectation value of a scalar field in the adjoint representation with the 
same form as Y: 



/2 



(9) 



-3/ 



The unbroken generators are those which commute with i.e. precisely the generators of SU{3) x SU{2) x U{1) 
above. 
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a) 





c) 




FIG. 1 Interference between the tree-level (a) and one-loop (b) diagrams with complex Yukawa couplings can provide the 
requisite source of CP violation for GUT baryogenesis. In viable models, to avoid the unwanted cancellations, one must often 
assume that the two scalars are different or go to higher loops (c) (Barr. Scgrc and Wcldon. .l979. : ...Kolb and Turner.. .1990.1 . 



The vector bosons which correspond to the broken generators gain mass of order gv. We will refer to the corre- 
sponding gauge bosons as X; they are associated with generators which don't commute with < <i> >, such as: 



/O 1 0\ 

' ' 



1 
Vo o/ 



(10) 



They carry color and electroweak quantum numbers and mediate processes which violate baryon number. For example, 
examining the structure of01 one sees that there is a coupling of the X bosons to a c? quark and an electron. Similarly, 
there is a coupling of the X boson to a quark doublet and a positron. Note that there is no way to assign baryon and 
lepton number to the X boson so that it is conserved by these couplings. 

In the GUT picture of baryogenesis, it is usually assumed that at temperatures well above the GUT scale, the 
universe was in thermal equilibrium. As the temperature drops below the mass of the X bosons, the reactions which 
produce the X bosons are not sufficiently rapid to maintain equilibrium. The decays of the X bosons violate baryon 
number; they also violate CP. So all three conditions are readily met: baryon number violation, CP violation, and 
departure from equilibrium. 

To understand in a bit more detail how the asymmetry can come about, note that CPT requires that the total 
decay rate of X is the same as that of its antiparticle X . But it does not require equality of the decays to particular 
final states (partial widths). So starting with equal numbers of X and X particles, there can be a slight asymmetry 
between the processes 

X^dL;X^Qu (11) 

and 

X^dl-X-^Qu. (12) 

The tree graphs for these processes are necessarily equal; any CP violating phase simply cancels out when we take 
the absolute square of the amplitude. This is not true in higher order, where additional phases associated with real 
intermediate states can appear. Actually computing the baryon asymmetry requires a detailed analysis, of a kind we 
will encounter later when we consider leptogenesis. 

There are reasons to believe, however, that GUT baryogenesis is not the origin of the observed baryon asymmetry. 
Perhaps the most compelling of these has to do with inflation. Assuming that there was a period of inflation, any 
pre-existing baryon number was greatly diluted. So in order that one produce baryons through X boson decay, it is 
necessary that the reheating temperature after inflation be at least comparable to the X boson mass. But as we have 
explained, a reheating temperature greater than 10^ GeV leads to cosmological difficulties, especially overproduction 
of gravitinos. 



D. Electroweak Baryon Number Violation 

Earlier, we stated that the renormalizable interactions of the Standard Model preserve baryon number. This 
statement is valid classically, but it is not qui te true of the c^uantum theory. There are, as we will see in this section, 
very tiny effects which violate baryon number l|'t Hooftlll976() . These effects are tiny because they are due to quantum 
mechanical tunneling, and are suppressed by a barrier penetration factor. At high temperatures, there is no such 
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suppression, so baryon number violation is a rapid process, which can come to thermal equilibrium. This has at 
least two possible implications. First, it is conceivable that these "sphaleron" processes can themselves be responsible 
for g enerating a baryon asymmetry. This is called electroweak baryogenesis IjKuzmin. Rubakov. and Shar)oshnikovi 
[l983). Second, as we will see, sphaleron processes can pr ocess an existing lepton number, producing a net lepton and 
baryon number. This is the process called leptogenesis (fFukugita and YanagidaL[l986j) . 

In this section, we summarize the main arguments that the electroweak interactions violate baryon number at high 
temperature. In the next section, we explain why the electroweak interactions might produce a small baryon excess, 
and why this excess cannot be large enough to account for the observed asymmetry. 

One of the great successes of the Standard Model is that it explains the observed conservation laws. In particular, 
there are no operators of dimension four or less consistent with the gauge symmetries which violate baryon number 
or the separate lepton numbers. The leading operators which can violate can baryon number are of dimension six, 
and thus suppressed by ©(-jp-). The leading operators which violate the separate lepton numbers are of dimension 
five, and thus suppressed by one power of 1/M. In each case, M should be thought of as the scale associated with 
some very high energy physics which violates baryon or lepton number. It cannot be determined except through 
measurement or by specifying a more microscopic theory. 

However, it is not quite true that the standard model preserves all of these symmetries. There are tiny effects, of 
order 

which violate them. These effects are related to the fact that the separate baryon number and lepton number currents 
are "anomalous." When one quantizes the theory carefully, one finds that the baryon number current, j^, is not 
exactly conserved but rather satisfies: 

Here F^^^, are the SU{2) field strengths, and we have introduced matrix-valued fields in the last expression, 

a 

and similarly for other fields, and the dual of F, F, is defined by: 

In electromagnetism, FF — 2E ■ B. 

The same anomaly ()14|l appears in the lepton number current as well, i.e., 

3-3 

^A'.^L ^ mrr^^f^"^'^'' ^ ^TrF^^F^^. (17) 

However, the difference of the two, jg — j^, is anomaly- free and is an exactly conserved quantity in the Standard 
Model (as well as SU(5) and SO(IO) grand unified theories). 

One might think that such a violation of current conservation would lead to dramatic violations of the symmetry. 
But the problem is more subtle. The right hand side of the anomaly equation is itself a total divergence: 

TrF^.F^, = dpK'' (18) 

where 



= e>'''''"tr[F,pA, + -A,ApA„] (19) 



(the reader can quickly check this for a U{1) gauge theory like electromagnetism). In view of this, 

J-fs-^K'^ (20) 

is conserved. In perturbation theory (i.e. in Feynman diagrams), falls to zero rapidly (typically like l/r^) at cx), 
and so its integral is zero. This fact insures that baryon number is conserved. 
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FIG. 2 Schematic Yang-Mills vacuum structure. At zero temperature, the instanton transitions between vacua with different 
Chern-Simons numbers are suppressed. At finite temperature, these transitions can proceed via sphalerons. 



In abelian gauge theories, this is the end of the story. In non-abelian theories, however, there are non-perturbative 
field configurations which contribute to the right hand side. These lead to violations of baryon number and the 
separate lepton numbers proportional to e~"^. Th ese configuration s are called instantons. We will not discuss them 
in detail here; a pedagogical treatment is given bv IColemanI l)l989f) . They correspond to calculation of a tunneling 
amplitude. To understand what the tunneling process is, one must consider more carefully the ground state of the 
field theory. Classically, the ground states are field configurations for which the energy vanishes. The trivial solution 
of this condition is A = 0, where A is the vector potential. More generally, one can consider A which is a "pure 
gauge," 

A=-g-^Vg, (21) 
I 

where 5 is a gauge transformation matrix. In an abelian (U(l)) gauge theory, fixing the gauge eliminates all but the 
trivial solution, A — Q} This is not the case for non-abelian gauge theories. There is a class of gauge transformations, 
labeled by a discrete index n, which do not tend to unity as \x\ —> cx3, which must be considered to be distinct states. 
These have the form: 

g^{x) = e»"/(^)*-^/2 (22) 

where f{x) —^2'Kasx—^ 00, and f{x) — > as x — > 0. 

So the ground states of the gauge theory are labeled by an integer n. Now if we evaluate the integral of the current 
K°, we obtain a quantity known as the Chern-Simons number: 



^ I d'xK" = d^xe,,kTrig-'d,gg-'d,gg~'dkg). 



(23) 



For g = gn, n^^ = n. The reader can also check that for g' = gn{x)h{x), where his a gauge transformation which tends 
to unity at infinity (a so-called "small gauge transformation"), this quantity is unchanged, n^^, the "Chern-Simons 
number," is topological in this sense (for A's which are not "pure gauge," n^^ is in no sense quantized). 

Schematically, we can thus think of the vacuum structure of a Yang-Mills theory as indicated in Fig. [21 We have, at 
weak coupling, an infinite set of states, labeled by integers, and separated by barriers from one another. In tunneling 
processes which change the Chern-Simons number, because of the anomaly, the baryon and lepton numbers will 
change. The exponential suppression found in the instanton calculation is typical of tunneling processes, and in fact 
the instanton calculation which leads to the result for the amplitude is nothing but a field-theoretic WKB calculation. 

At zero temperature, the decay amplitude is suppressed, not only by , but by factors of Yukawa couplings. 
The probability that a single proton has decayed through this process in the history of the universe is infinitesimal. 
But this picture suggests that, at finite temperature, the rate should be larger. One can determine the height of the 



More precisely, this is true in axial gauge. In the gauge Ao = 0, it is necessary to sum over all time-independent transformations to 
construct a state which obeys Gauss's law. 
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barrier separating configurations of different n^^ by looking for the field configuration which corresponds to sitting 
on top of t he barrier. Th is is a solution of the static equations of motion with finite energy. It is known as a 
"sphaleron" l|Mantoni Il98^ . When one studies the small fluctuations about this solution, one finds that there is a 
single negative mode, corresponding to the possibility of rolling down hill into one or the other well. The sphaleron 
energy is of order 

E,p = -^M„. (24) 

This can be seen by scaling arguments on the classical equations; determining the coefficient c requires a more 
detailed analysis. The rate for thermal fiuctuations to cross the barrier per unit time per unit volume should be of 
order the Boltzmann factor for this configuration, times a suitable orefactor l, Arnold and McLerran. .lOSR: ..Dine et.al . 
uzmin. Rubakov. and Shaposhnikovlll985|) . 

Tsp = T^e"^-/^. (25) 

Note that the rate becomes large as the temperature approaches the W boson mass. In fact, at some temperature 
the weak interactions undergo a phase transition to a phase in which the W boson mass vanishes. At this point, the 
computation of the transition rate is a difficult problem - there is no small parameter - but general scaling arguments 
show that the transition rate is of the form^ : 

Tbv = al,T\ (26) 

Returning to our original expression for the anomaly, we see that while the separate baryon and lepton numbers 
are violated in these processes, the combination B — L is conserved. This result leads to three observations: 

1. If in the early universe, one creates baryon and lepton number, but no net B — L, B and L will subsequently 
be lost through sphaleron processes. 

2. If one creates a net B — L (e.g. creates a lepton number) the sphaleron process will leave both baryon and 
lepton numbers comparable to the original B — L. This realization is crucial to the idea of leptogenesis, to be 
discussed in more detail below. 

3. The standard model satisfies, by itself, all of the conditions for baryogenesis. 



E. Electroweak baryogenesis 

A s we will see, while the Standard Model satisfies all of the conditions for baryogene- 
sis l)Kuzmin. Rubakov. and Shaposhnikovl Il985l) . nothing like the required baryon number can be produced. 
It is natural to ask whether extensions of the Standard Model, such as theories with complicated Higgs, or the 
Minimal Supersymmetric (extension of the) Standard Model, can generate an asymmetry, using the sphaleron process 
discussed in the previous section. We will refer to such a possibility more generally as "Electroweak Baryogenesis." 



1. Electroweak baryogenesis in the Standard Model 

How might baryons be produced in the Standard Model? From our discussion, it is clear that the first and second 
of Sakharov's conditions are satisfied. What about the need for a departure from equilibrium? 

Above we alluded to the fact that in the electroweak theory, there is a phase transition to a phase with massless 
gauge bosons. It turns out that, for a sufficiently light Higgs, this transition is first order. At zero temperature, in 
the simplest version of the Standard Model with a single Higgs field, 0, the Higgs potential is given by 

Vm=-ii^m^ + ^\^\ (27) 

The potential has a minimum at $ = -^Vo: breaking the gauge symmetry and giving mass to the gauge bosons by 
the Higgs mechanism. 



^ more detailed considerations alter slightly the parametric form of the rate iArnold. Son, and Yaffel 119971) 
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FIG. 3 First and second order phase transitions. 

What about finite temperatures? By analogy with the phase transition in the Landau- Ginsburg model of super- 
conductivity, one might expect that the value of < $ > will change as the temperature increases. To determine the 
value of $, one must compute the free energy as a function of $.The leading temperature-dependent corrections are 
obtained by simply noting that the masses of the various fields in the theory - the W and Z bosons and the Higgs 
field, in particular, depend on <i>. So the contributions of each species to the free energy are ^-dependent: 



where (3 = l/T, T is the temperature, the sum is over all particle species (physical helicity states), and the plus sign 
is for bosons, the minus for fermions. In the Standard Model, for temperature T ~ 10^ GeV, one can treat all the 
quarks are massless, except for the top quark. The effective potential H28() then depends on the top quark mass, , 
the vector boson masses, and , and on the Higgs mass, . Performing the integral in the equation yields 



The parameters To, D and E are given in terms of the gauge boson masses and the gauge couplings below. For the 
moment, though, it is useful to note certain features of this expression. E turns out to be a rather small, dimensionless 
number, of order 10~^. If we ignore the (p^ term, we have a second order transition, at temperature To, between a 
phase with and a phase with </> = 0. Because the W and Z masses are proportional to (f>, this is a transition 

between a state with massive and massless gauge bosons. 

Because of the (j)^ term in the potential, the phase transition is potentially at least weakly first order. This is 
indicated in Fig. Here one sees the appearance of a second, distinct, minimum at a critical temperature. A first 
order transition is not, in general, an adiabatic process. As we lower the temperature to the transition tempera- 
ture, the transition proceeds by the formation of bubbles; inside the bubble the system is in the true equilibrium 
state (the state which minimizes the free energy) while outside it tends to the original state. These bubbles form 
through thermal flucutations at different points in the system, and grown until they collide, completing the phase 
transition. The moving bubble walls are regions where the Higgs fields are changing, and all of Sakharov's conditions 
are satisfied. It has been shown that various non-equilibrium processes n ear th e wall can produce baryon and lepton 
numbers (Cohen. Kaolan and Nelson, 1993; Rubakov and Shaposhnikovl [l996() . 

Describing these processes would take us far afield. Even without going through these details, however, one point is 
crucial: after the bubble has passed any given region, the baryon violating processes should shut off. If these processes 
continue, they wash out the baryon asymmetry produced during the phase transition. 

Avoiding the washing out of the asymmetry requires that after the phase transition, the sphaleron rate should be 
small compared to the expansion rate of the universe. According to eq. (|25|) . this requires that after the transition the 
sphaleron energy, which is proportional to (temperature-dependent) M^-boson mass, M^, be large compared to the 
temperature. This, in turn, means that the Higgs expectation value must be large immediately after the transition. 
Using eq. (|29fl or more refined calculations to higher orders, one can relate the change in the Higgs expectation value 




(28) 



y ($, T) = T)(T2 - T2)$2 _ ET<i>^ + + .... 



(29) 
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to the Higgs mass at zero temperature. It turns out that the current lower limit on the Higgs boson mass rules out any 
possibility of a large enough Higgs expectation value immediately after the phase transition, at least in the minimal 
model with a single Higgs doublet. 

The shape of the effective potential V{(l), T) near the critical temperature determines whether the phase transition 
is first order, which is a necessary condition for electroweak baryogenesis to work. E g. (I29II represents the lowest 
order term in perturbation theory; higher order terms hav e been computed as well l|Arnold and EspinosaL Il993t 
iBagnasco an d Dinel Il993t iFarakos et alV Il994t iLainel Il99^ . However, these calculations are not reliable, because 
of infrared divergences which arise in the perturbation expansion. These arise because the Higgs field is nearly 
massles s at the transition. Numerical simulations ar e required, o f ten combined with a clever use of perturbation 
theor y l|Fa rakos e t a/.l. 119951). Numerical simulation s llCsikor et al\. Il999t iGurtler et al\. Il997t iKaiantie et alV [1996', 
11993.11999; Karsch et a^J . ll997t [Rummukaine n~ei aZj . lT998[) have shown that, for the Higgs mass above 80 GeV (which 
it must be, to satisfy the present experimental constraints), the sharp phase transition associated with the low mass 
Higgs turns into a smooth crossover. 

However, even for an unrealistically light Higgs, the actual production of baryon asymmetry in the minimal Stan- 
dard Model would be highly suppresse d. This is because Standard Model CP violation must involve all three gener- 
ations ijKobavashi and Maskawaril973(l . The lowest order diagram that involves three generations of fermions with 
proper chi ralities and contributes to CP violating processes relevant to baryogenesis is suppressed by 12 Yukawa 
couplings ijShaposhnikovl Il986l 1198'^ . Hence, the CKM CP violation contributes a factor of 10^^" to the amount of 
baryon asymmetry that could arise in the Standard Model. 

Clearly, one must look beyond the Standard Model for the origin of baryon asymmetry of the universe. One of the 
best motivated candidates for new physics is supersymmetry. 

Before closing this section, for completeness, we give the values of the parameters To, B, D and E in can. 1291 

To = ^(m' - ^Bvl) = -^{ml - %Bvl) , (30) 
while the parameters -B, D and E are given by: 

^ = 70^(2^^ + ^^l-4m4), D^^{2M^ + ml + 2m]), £; ^ J^(2M3, + m|) ^ 10"^. (31) 



2. Supersymmetry, a short introduction 



In this section we provide a brief introduction to supersymmetry. Much more detail can be found, e.g., in IjDind . 
and in many texts. 

There are many hints that supersymmetry, a hypothetical symmetry between fermions and bosons, might play some 
role in nature. For example, supersymmetry seems to be an essential part of superstring theory, the only consistent 
theory of quantum gravity which we know. If supersymmetry is a symmetry of the laws of nature, however, it must 
be badly broken; otherwise we would have seen, for example, scalar electrons ( "selectrons" ) and fermionic photons 
( "photinos" ) . It has been widely conjectured that supersymmetry might be discovered by accelerators capable of 
exploring the TeV energy range. There are several reasons for this. The most compelling is the "hierarchy problem." 
This is, at its most simple level, the puzzle of the wide disparity of energies between the Planck scale (or perhaps 
the unification scale) and the weak scale - roughly 17 orders of magnitude. While one might take this as simply a 
puzzling fact, within quantum theory, the question is made sharper by the fact that scalar masses (particularly the 
Higgs mass) receive very divergent quantum corrections. A typical expression for the quantum corrections to a scalar 
mass is: 

Sm^^ - [ d^k^. (32) 
n J 

This integral diverges quadratically for large momentum (fc). Presumably, the integral is cut off by some unknown 
physics. If the energy scale of this physics is A, then the corrections to the Higgs mass are much larger than the 
scale of weak interactions unless A ^ TeV. While various cutoffs have been proposed, one of the most compelling 
suggestions is that the cutoff is the scale of supersymmetry breaking. In this case, the scale must be about 1000 GeV. 
If this hypothesis is correct, the Large Hadron Collider under construction at CERN should discover an array of new 
particles and interactions (it is possible that supersymmetry could be discovered at the Tevatron beforehand). 

The supersymmetry generators, Qa, are fermionic operators. Acting on bosons they produce fermions degenerate 
in energy; similarly, acting on fermions, they produce degenerate bosons. Their algebra involves the total energy and 
momentum. 



(33) 
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Neglecting gravity, supcrsymmctry is a global symmetry. Bccaiisc of the structure of the algebra, the symmetry is 
broken if and only if the energy of the ground state is non-zero. If the symmetry is unbroken, for every boson there 
is a degenerate fermion, and conversely. 

If we neglect gravity, there are two types of supermultiplets which may describe light fields. These are the chiral 
multiplets, containing a complex scalar and a Weyl (two-component) fermion 

= V'i), (34) 

and the vector multiplets, containing a gauge boson and a Weyl fermion: 

y« = (^«,A«). (35) 

In global supersymmetry, the lagrangian is specified by the gauge symmetry and an analytic (more precisely holo- 
morphic) function of the scalar fields, W{^i), known as the superpotential. For renormalizable theories, 



W{4>i) = -■mij(t>i4>j + \ijk<l>i(t>j4>k- (36) 



The lagrangian then includes the following: 

1. The usual covariant kinetic terms for all of the fields, for example 

A»£»A-, i,Dr, ^i^;.', \D^<i>?. (37) 

2. Yukawa couplings with gauge strength: 

\/25"A"(?!)*T«V + c.c. (38) 

3. Mass terms and Yukawa couplings from W: 

1 d'^W 



2, 



-i>i^j (39) 



4. A scalar potential: 



3 

= rriijipi'ipj + -Xijk4>2i^jipk- (40) 



It is convenient to define two types of auxiliary fields, the "F" and "D" fields: 

dW 

d(j)i 



^' = ^ i5<^ = 9"^<^*T>,. (42) 



In terms of these, the potential is simply 



F= + (43) 



and, at the classical level, supersymmetry is unbroken if and only if all of the D and F fields vanish at the minimum 
of the potential. 

It is useful to consider some examples. Take first a model with a single chiral field, (j), and superpotential 

W = ^m<tP-. (44) 
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In this case, the potential is 

2 

= m2|(/,|2. (45) 

On the other hand, the fermion mass is just m, from eqn. H4U|) . so the model describes two bosonic and two fermionic 
degrees of freedom, degenerate in mass. 

A more interesting model is a supersymmetric version of the standard model, known as the Minimal Supersymmetric 
Standard Model (MSSM). The gauge group is SU{3) x SU{2) x U{1), and there is one vector multiplet for each gauge 
generator. In addition, for each of the usual fermions of the standard model, one has a chiral field with the same 
quantum numbers: 

Qa = (3,2)i/3 S, = (3,2)_4/3 4 = (3,2)2/3 -(l,2)_i e = (1,1)2. (46) 

Here a is a generation index, a = 1, ... ,3. In addition, there are two Higgs fields (two in order to cancel anomalies 
and to be able to give mass to all quarks and leptons) 

i/„ = (1,2)1 i/d = (l,2)_i (47) 

In this notation, (x, 'iI')y shows the transformation properties of the fields under the action of the gauge group: x ^^'^ 
ip refer to the representations of the SU(3) and SU(2) respectively, while Y denotes the hypercharge. 
The superpotential of the model is a generalization of the Yukawa couplings of the standard model: 

W = T'^'^QaUbH^ + l^QadaHd + TaLaEaHa + flHuHd, (48) 

where we have used our freedom to make field redefinitions to make the (somewhat unconventional, but later conve- 
nient) choice that the c?-quark and the lepton Yukawa couplings are diagonal. If one supposes that the Higgs fields 
Hu and have expectation values, this gives, through equation (|39|l . masses for the quarks and leptons just as in 
the Standard Model. If supersymmetry is unbroken, their scalar partners obtain identical masses. The last term is a 
supersymmetric mass term for the Higgs fields. Supersymmetry breaking, essential to obtain a realistic model, will 
be discussed momentarily. 

The gauge symmetries actually permit many more couplings than those written in eqn. iflSI) . Couplings such as 
HL, udd, and others would violate baryon or lepton number if they appeared. Because these are dimension four, they 
are unsuppressed (unless they have extremely tiny dimensionless coefficients) . They can be forbidden by a symmetry, 
under which ordinary fields are even (quarks, leptons, and Higgs bosons) while there supersymmetric partners are 
odd. This symmetry is called R-parity. 

By itself, this model is not realistic, since supersymmetry is unbroken and all ordinary fields (quarks, leptons, gauge 
bosons, higgs) are degenerate with their superpartners (squarks, sleptons, gauginos). The simplest solution to this is 
just to add "soft breaking terms" which explicitly break the supersymmetry. Because they are soft, they don't spoil 
the good features of these theories. These soft terms include mass terms for the squarks and sleptons, Majorana mass 
terms for the gauginos, and cubic couplings of the scalar fields, 

ml\<f>"^ P + mxXX + mAA,,kC^'(jy>(j)'' (49) 

In the minimal supersymmetric standard model, there are 105 such couplings (counting real parameters). We will 
think of all of these mass parameters as being of order 7713/2 ~ "niz. These parameters are highly constrained, both 
by low energy physics (particularly by the suppression of flavor-changing processes in weak interactions) and direct 
searches at LEP and the Tevatron. Theoretical approaches to understanding these soft breakings can be divided 
broadly into two classes. Both assume that some dynamics gives rise to spontaneous breakdown of supersymmetry. 
In "gravity mediation," very high energy physics is responsible for generating the soft terms; in gauge mediated 
models, lower energy, gauge interactions communicate supersymmetry breaking to ordinary fields. 

R parity, if present, implies that the lightest of the new particles, called the "LSP" implied by supersymmetry is 
stable. Typically this is the partner of a neutral gauge or Higgs boson. One can calculate the abundance of these 
particles (neutralinos) as a function of the various supersymmetry breaking parameters. The assumption that the 
supersymmetry-breaking masses are hundreds of GeV leads automatically to a neutralino density of order the dark 
matter density of the universe, and this particle is a leading candidate for the dark matter. 

3. Baryogenesis in the MSSM and the NMSSM 

Supersymmetric extens i ons of the Standard Model contain new sources of CP violation ijDine et.al. Il99lt 
iDine. Huet. and Singletonl Il992t iHuet and NelsonL Il996j) and an enlarged set of parameters which allow a greater 
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possibility of a first-order trans ition ("Bodckcr et aL'.'l QQTtlCarena. Quiros. and Waeneil 11998*: 'Cli ne and MooreLll998t 
lEspinosa. Quiros. and Zwirneil 11993; .Espinosa. .19961) . So it would seem possible that electroweak baryogenesis could 
operate effectively in these theories. 

The new sources of CP violation may come, for example, from the chargino mass matrix: 

= sj), ( ^ ) ( |; )^ + (50) 

where w and h are the superpartners of W-boson and the charged Higgs, and TO2 and /i are some mass parameters. 

Other possible sources include phases in scalar masses (by field redefinitions, some of these can be shifted from 
fermion to scalar mass terms). We will focus, however, on the terms in eq. (I5U|) . 

As long as TO2 and u are complex, spatially varying phases in the bubbl e wall provide a source of (spontaneous) 
CP violation ICohen. Kaplan, and Nelstml Il99lt lO Il974t IWeinberd. I1976D . However, in light of the constraints on 
Higgs and s uperpartner masses, the present w indow f or electroweak baryogenesis in the M SSM is very narrow if it 
exists at all (jCarena. Quiros. Seco. and Wagn er. 2003 lCline. Jovce and KainulainenL Il998() . As we discussed above, 
the lighter the Higgs the easier it is to avoid the wash-out of baryon asymmetry produced in the phase transition. A 
light right-handed stop allows for a first-order phase transition even for Higgs as heavy as 115 GeV, which is barely 
consistent with the current bounds. The predictions of the light Higgs and the light stop will soon be tested in 
experiment. 

However, even for the most optimistic choice of parameters, it is difficult to obtain a baryon asymmetry as large as 
the observed value quoted in eq. Several parameters must be adjusted to maximize the baryon asymmetry. In 
particular, one must assume that the wall is very thin and choose the "optimal" bubble wall velocity v^, ~ 0.02. The 
origin of these difficulties lies, once again, in the strength of the electroweak phase transition. In the MSSM, the phase 
transition can be enhanced if the right-h anded stop (the scalar partner of the top quark) is assumed t o be very light , 
while the left-handed stop is very heavy ijCarena. Quiros. and WagneRll998() . Then two- loop effects ijBodeker et all 
I1997I: lEspinosal Il996|) change the scalar pote ntial sufficiently to allow for a first-order phase tra nsition; the lattice 
simulations support this perturbative result ijCsikor et all l2000t II aine and Rummukainerj . [199^. However, severe 
constraints arise from the experiment al bounds on the chargino mass, as well as the chargino contribution to the 
electric dipole moment of the neutron ijChang. Chang and KeunEll2002t|Pilaftsisll2002(l . 

Different calculations of the baryon asymmetry in the MSSM yie ld somewhat different re- 
sults (Ca rena. Quiros, Seco, and Wagner, 2002; Cline, Joyce and Kainulainen, .1998(1 . as can be seen from Fig. 0] 
According to lCarena. Quiros. Seco." and Wagner! |200^, it is possible to produce enough baryons if the Higgs boson 
and the right-handed stop are both very light, near the present experimental limits. In any case, electroweak baryo- 
genes is in the MSSM is on the verge of being confirmed or ruled out by improving experimental constraints fClinel 
120021) . 

The strength of the phase transition can be further enhanced by adding a singlet Higgs to the model. 
In th e next-to- minimal supersymmetr i c model (NMSSM), the phase transition can be more strongly first- 
order l|Davies. Froggatt. and Moorhouselll996tlHuber et a^J . l200 lllKainu laincn et al'.. '2OOII. The singlet also provides 
additional sources of CP violation which increase the baryon asymmetry (,Huber and Schmidt . .2001.1 . 

4. Non-thermal electroweak baryogenesis at preheating 

In light of these difficulties, various proposals have been put forth to obtain a viable picture of electroweak baryo- 
genesis. These typically involve more drastic departures from thermal equilibrium than the weakly first order phase 
transitions described above. The more extreme proposals suppose that inflation occurred at the electroweak scale 
and kicked the universe out of equilibrium, setting the stage for baryogenesis. It i s generally believed that the nat- 
ural scale for inflation is m uch higher than 10^ GeV. Although models with we ak l|German. Ross, and Sarkail l200lt 
iRandall and Thomaslll99l or intermediate jRandall. Soliacic. and Guthl . ll996(l scale inflation have been constructed, 
a lower scale of inflation is generally difficult to reconcile with the observed density perturbations (Sp/p) ~ 10"^. 
As a rule, the smaller the scale of inflation, the flatter the infiaton potential must be to produce the same density 
fluctuations. A weak-scale inflation would require the infiaton potential to be extremely flat, perhaps flatter than can 
plausibly be obtained in any physical theory. 

Qf course, one does not have to assume that the same inflation is responsible for (5p/ p) and for baryogenesis. Qne 
could imagine that the universe has undergone more than one inflationary period. The primary inflation at a high 
scale could be responsible for the flatness of the universe and for the observed density perturbations. A secondary 
inflation at the weak scale need not produce an enormous expansion of the universe, and could create fertile soil 
for baryogenesis. Qne can debate the plausibility of invoking a second stage of infiation just for this purpose. In 
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FIG. 4 Results of two independent calcul ations of baryon asymmetry in the M SSM. On the left, contours of constant baryon 
asymmetry in units 10~^° calculated by ICline. Joyce and KaimilainenI (^^E) for the bubble wall velocity «w = 0.03 and 
tan/3 ~ 3. Mass units are GeV. Shaded regions are excluded by the LEP2 limit on the chargino mass, m^± > 104 GeV. To 
maximize the baryon asymmetry, one assumes that the bubble wall is very narrow, ^ 6/T, and its velocity is — 0.03. 
The plot on the right represents the results of ICarena. Quiros. Seco. and Wagneil fcOOST) for tan/3 = 10 and a mgiximal CP 
violating phase of the fi parameter. The other parameters in both calculations have been chosen to maximize the resultin g 
baryon asymmetry. Also shown is the observed value of the baryon asymmetry reported by WMAP iBennett et all 120031) , 

„ _ f! 1 +0.3 
'/ — "-J^ -0.2- 

favor of such a possibility, it has been argued that a low-scale inflation might ameliorate the cosmoloKical modu li 
problem common to many supersymmetric theories (jGerman. Ross, and SarkaR l2nOH iRandall and Thomal llQOfJ) . 
Nevertheless, inflationary models at th e electroweak scale, which would help generate baryon asymmetry, usually 
suffer from naturalness problems ijLvthl . Il9991. which may be less severe in some cases ( (Cop eland et ai, 2001). 

What inhibits electroweak baryogenesis in the Standard Model is too much equilibrium and too little CP violation. 
Both of these problems might be rect ified if in flation is followed by re heating to a temperature just below the 
electroweak scale ( Garcia-Bclh do et all |l999; Kr auss and TroddenL Il999^■ Reheating, especially its vari ant dubbed 
preheating, involves a radical departure from thermal equilibrium ('Kof man. Linde and Starobinskvl Il996() . 

During inflation, all matter and radiation are inflated away. When inflation is over, the energy stored in the 
inflaton is converted to thermal plasma. There are several possibilities for this reheating process. One possibil- 
ity is that the inflaton may decay perturbatively into light particles, which eventually thermalize. However, in a 
class o f models, a parametric resonance may greatly enhance the production of particles in some specific energy 
bands ijKofman. Linde and Starobinskvlll996j) . This process, caused by coherent oscillations of the inflaton, is known 
as preheating. Alternatively, the motion of the condensate may become spatially inhomogenious on scales smaller 
than the horizon. This kind of tr ansition in the mot ion of the inflaton, called spinodal decomposition, may lead to a 
very rapid tachyonic preheating ijFelder et al\ . \200^ . 

All of these variants of reheating force the universe into a non-equilibrium state after the end of inflation and 
before thermalization takes place. This is, obviously, an opportune time for baryogenesis. The usual considerations of 
sphaleron transitions do not apply to a non-equilibrium system. But it turns out that baryon-number violating pro- 
cesses similar to sphaleron transitions do take place at preheating ijCornwall and Kusenkol 120001 : iGarcia-Bellido et all . 
^99), as well as tachyonic preheating ijGarcia-Bellido et aZ.L l2003t ISmit and TranbereLl200^ . This has been demon- 
strated by a combination of numerical and analytical a rguments. In ad dition, preheating allows the coherent motions 
of some condensates to serve as sources of CP violation (jCornwall et a/.l . [2001.) . Such sources are poorly constrained by 
experiment and could have significant impact on baryogenesis. It is conceivable, therefore, that the electroweak-scale 
inflation could facilitate generation of the baryon asymmetry. 

F. Leptogenesis 

Of the five scenarios for baryogenesis which we have listed in the introduction, we have discussed two which are 
connected to very high energy physics: Gravitational and GUT baryogenesis. We have given cosmological arguments 
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why they are not hkely. These arguments depend on assumptions which we cannot now rehably estabhsh, so it is 
yet possible that these mechanisms were operative. But if we tentatively accept these arguments we can significantly 
narrow our focus. Similarly we have seen that electroweak baryogenesis, while a beautiful idea, cannot be implemented 
in the Standard Model, and probably not in its minimal supersymmetric extension. So again, while we can not rule 
out the possibility that electroweak baryogenesis in some extension of the Standard Model is relevant, it is tempting, 
for the moment, to view this possibility as unlikely. Adopting this point of view leaves leptogenesis and Affleck-Dine 
baryogenesis as the two most promising possibilities. What is exciting about each of these is that, if they are operative, 
they have consequences for experiments which will be performed at accelerators over the next few years. 

While there is no experimental evidence for supersymmetry apart from the unification of couplings, in the l ast few 
years, evidence for neutrino masses has become more and more compelling ijGonzalez-Garcia and NiJ . |2002|) . This 
comes from several sources: the fact that the flux of solar neutrinos does not match theoretical expectations, in the 
absence of masses and mixings ( Super- Ki l200lt FSNOl 120021) : the apparent observation of neutrino oscil lations among 
atmospheric neutrinos (jSuper-K. .1998^) : and direct measurements of neutrino mixing llKamLANDl l2002). 

We will not review all of these phenomena here, but just mention that the atmospheric neutrino anomaly suggests 
oscillations between the second and third generation of neutrinos: 

Ami = 10"^ - lO^W (51) 

with mixing of order one, while the solar neutrino deficit suggests smaller masses (Am2 - 10-^eV^). There is other 
evidence for neutrino oscillation from accelerator experiments. The SNO experiment has recently provided persuasive 
evidence in support of the hypothesis of mixing (as opposed to modifications of the standard solar model). The 
results from Super-Kamiokande, SNO and KamLAND are in good agreement. There is also evidence of mixing from 
an experiment at Los Alamos (LSND). This result should be confirmed, or not, by the MiniBoone experiment at 
Fermilab. The mixing suggested by atmospheric neutrinos is currently being searched for directly by accelerators. 
The data so far supports the mixing interpretation, but is not yet decisive. 

The most economical explanation of these facts is that neutrinos have Majorana masses arising from lepton-number 
violating dimension five operators. (A Majorana mass is a mass for a two component fermion, which is permitted if 
the fermion carries no conserved charges.) We have stressed that the leading operators permitted by the symmetries of 
the Standard Model which violate lepton number are non-renormalizable operators of dimension five, i.e. suppressed 
by one power of some large mass. Explicitly, these have the form: 

Ci, = ^LHLH (52) 

2 

Replacing the Higgs field by its expectation value v gives a mass for the neutrino of order jj. If M = Mp, this 
mass is too small to account for either set of experimental results. So one expects that some lower scale is relevant. 
The "see-saw" mechanism provides a simple picture of how this scale might arise. One supposes that in addition 
to the neutrinos of the Standard Model, there are some SU{2) x [/(l)-singlet neutrinos, N. Nothing forbids these 
from obtaining a large mass. This could be of order M^^^, for example, or a bit smaller. These neutrinos could also 
couple to the left handed doublets i^^^ , just like right handed charged leptons. Assuming, for the moment, that these 
couplings are not particularly small, one would obtain a mass matrix, in the {N, v^} basis, of the form 

This matrix has an eigenvalue -jf^. The latter number is of the order needed to explain the neutrino anomaly for 
M ~ IG^'^ or so, i.e. not wildly different than the GUT scale and other scales which have been proposed for new 
physics. 

For leptogenesis l)Fukugita and Yanagidal|l986fl . what is important in this model is that the couplings of TV break 
lepton number. iV is a heavy particle; it can decay both to h + h' and h + for example. The partial widths to 
each of these final states need not be the same. CP violation can enter through phases in the Yukawa couplings 
and mass matrices of the iV's. At tree level, however, these phases will cancel out between decays to the various 
states and their (would be) CP conjugates, as in the case of GUTs we discussed earlier. So it is necessary to consider 
interference between tree and one loop diagrams with discontinuities, as in Fig.^ In ^ model with three TV's, there are 
CP-violating phases in the Yukawa couplings of the iV's to the light Higgs, The heaviest of the right handed neutrinos, 
say A^i, can decay to £ and a Higgs, or to £ and a Higgs. At tree level, as in the case of GUT baryogenesis, the rates 
for production of leptons and antileptons are equal, even though there are CP violating phases in the couplings. It is 
necessary, again, to look at quantum corrections, in which dynamical phases can appear in the amplitudes. At one 
loop, the decay amplitude for iV has a discontinuity associated with the fact that the intermediate Ni and N2 can 
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be on shell. So one obtains an asymmetry proportional to the imaginary parts of the Yukawa couplings of the iV's to 
the Higgs: 



i=2,3 ^ 1 



(55) 



where / is a function that represents radiative corrections. For example, in the Standard Model / = \/x[{x ~ 2)/ {x — 
1) + (a; + 1) ln(l + l/x)], while in the MSSM / = ^[2/{x - 1) + ln(l + l/x)]. Here we have allowed for the possibility 
of multiple Higgs fields, with H2 coupling to the leptons. The rough order of magnitude here is readily understood 
by simply counting loops factors. It need not be terribly small. 

Now, as the universe cools through temperatures of order the of masses of the iV's, they drop out of equilibrium, and 
their decays can lead to an excess of neutrinos over antineutrinos. Detailed predictions can be obtained by integrating 
a suitable set of Bo ltzman equations. Alternatively, th ese particles can be produced out of equilibrium, at prehearing 
following inflation IjGarcia- BeUido and Moralesl I200H1 . 

These decays, then, produce a net lepton number, but not baryon number (and hence a net B — L). The resulting 
lepton number will be further processed by sphaleron interactions, yielding a net lepton and baryon number (recall that 
sphaleron interactions preserve B — L, but violate B and L separately ). One can determine the resulting asymmetry 
by an elementary thermodynamics exercise l|Harvev and Turneil llQQCl) . One introduces chemical potentials for each 
neutrino, quark and charged lepton species. One then considers the various interactions between the species at 
equilibrium. For any allowed chemical reaction, the sum of the chemical potentials on each side of the reaction must 
be equal. For neutrinos, the relations come from: 



1. the sphaleron interactions themselves: 



2. a similar relation for QCD sphalerons: 



^(3^,, + = (56) 



^(2^q, - - /irfj = 0. (57) 

i 

3. vanishing of the total hypercharge of the universe: 

(^9' ~ + ^rf. ~ + Me.) + = (58) 

i 

4. the quark and lepton Yukawa couplings give relations: 

A*9i - M0 - Mdj = 0- Mg. - M0 - = 0, - M0 ~ Me^ = 0. (59) 

The number of equations here is the same as the number of unknowns. Combining these, one can solve for the 
chemical potentials in terms of the lepton chemical potential, and finally in terms of the initial B — L. With N 
generations, 

B^ — (B-L). (60) 

22Ar+13^ ' ^ ' 

Reasonable values of the neutrino parameters give asymmetries of the order we seek to explain. Note sources of 
small numbers: 

1. The phases in the couplings 

2. The loop factor 

3. The small density of the N particles when they drop out of equilibrium. Parametrically, one has, e.g., for 
production, 

r ^ e^-^t/T)g^T (61) 

which is much less than H ^ T'^/Mp once the density is suppressed by T/Mp, i.e. of order 10~^ for a 10"'^^ GeV 
particle. 
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It is interesting to ask: assuming that these processes are the source of the observed asymmetry, how many 
parameters which enter into the computation can be measured, i.e. can we relate the observed number to microphysics. 
It is hkely that, over time, many of the parameters of the li ght neutrino mass matrices, including possible CP-violating 
effects, will be measured llGonzalez-Garcia and Nill2002|) . But while these measurements determine some of the Ni 
couplings and masses, they are not, in general, enough. In order to give a precise calculation, analogous to the 
calculations o f nucleosynthesis, of the baryon number density, one needs additional information about the masses of 
the fields iPa, sco]i p.t oil 12003^. One either requires some other (currently unforseen) experimental access to this 
higher scale physics, or a compelling theory of neutrino mass in which symmetries, perhaps, reduce the number of 
parameters. 

G. Baryogenesis through Coherent Scalar Fields 

We have seen that supersymmetry introduces new possibilities for electroweak baryogenesis. But the most striking 
feature of supersymnietric models, from the point of view of baryogenesis, is the appearance of scalar fields carrying 
baryon and lepton number. These scalars offer the possibility of coherent production of baryons. In the limit that 
supersymmetry is unbroken, many of these scalars have flat or nearly flat potentials. They are thus easily displaced 
from their minima in the highly energetic environment of the early universe. We will often refer to such configurations 
as "excited." Simple processes can produce substantial amounts of baryon number. This coherent production of 
baryons, known as Affleck-Dine baryogenesis, is the focus of the rest of this review. 

III. AFFLECK-DINE BARYOGENESIS 

A. Arguments for Coherent Production of the Baryon Number 

In the previous section, we have reviewed several proposals for generating the baryon number. None can be firmly 
ruled out, however all but two seem unlikely: leptogenesis and Affleck-Dine baryogenesis. While the discovery of 
neutrino mass gives support to the possibility of leptogenesis, there are a number of reasons to consider coherent 
production: 

• The Standard Model alone cannot explain the baryon asymmetry of the universe, the main obstacle being the 
heaviness of the Higgs. One needs new physics for baryogenesis. The requisite new physics may reside at a 
very high scale or at a lower scale. An increasing body of evidence implies that inflation probably took place 
in the early universe. Hence, baryogenesis must have happened at or after reheating. To avoid overproducing 
weakly interacting light particles, for example gravitinos and other new states predicted in theory, one would 
like the reheat temperature not to exceed 10^ GeV. This poses a problem for GUT baryogenesis. This also limits 
possibilities for the leptogenesis. Affleck-Dine baryogenesis, on the other hand, is consistent with low energy 
and temperature scales required by inflation. 

• Supersymmetry is widely regarded as a plausible, elegant, and natural candidate for physics beyond the Standard 
Model. Of the two simple scenarios for baryogenesis in the MSSM, the electroweak scenario is on the verge of 
being ruled out by accelerator constraints on supersymmetric particles, in sharp contrast with the AD scenario. 

• The remaining low-reheat SUSY scenario, Affleck-Dine baryogenesis, can naturally reproduce the observed 
baryon asymmetry of the universe. The formation of an AD condensate can occur quite generically in cosmo- 
logical models. 

• The Affleck-Dine scenario potentially can give rise simultaneously to the ordinary matter and the dark matter 
in the universe. This can explain why the amounts of luminous and dark matter are surprisingly close to each 
other, within one order of magnitude. If the two entities formed in completely unrelated processes (for example, 
the baryon asymmetry from leptogenesis, while the dark matter from freeze-out of ncutralinos), the observed 
relation JIdark ^ f^mattcr is fortuitous.'^ 



^ An additional ad hoc symmetry can also help relate the amounts of ordinary matter and dark matter ^KaDlanLll992^ . 
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• Many particle physics models lead to significant production of entropy at relatively late 
times l|Cohen. Kaplan and NelsonL Il993() . This dilutes whatever baryon number existed previously. Co- 
herent production can be extremely efficient, and in many models, it is precisely this late dilution which yields 
the small baryon density observed today. 

In the rest of this section, we discuss AfBeck-Dine baryogenesis in some detail. 



B. Baryogenesis Through a Coherent Scalar Field: 

In supersymmetric theories, the ordinary quarks and leptons are accompanied by scalar fields. These scalar fields 
carry baryon and lepton number. A coherent field, i.e., a large classical value of such a field, can in principle carry a 
large amount of baryon number. As we will see, it is quite plausible that such fields were excited in the early universe. 

To understand the basics of the mechanism, consider first a model with a single complex scalar field. Take the 
lagrangian to be 

£= |a^0|2-m2|0|2 (62) 

This lagrangian has a symmetry, (j) — > e*"^, and a corresponding conserved current, which we will refer to as baryon 
number: 



It also possesses a "CP" symmetry: 



j>^^i{(t)*d^'(j)-(t)d^'(t)*). (63) 
(^^(/)*. (64) 



With supersymmetry in mind, we will think of m as of order M^y. 

If we focus on the behavior of spatially constant fields, (f>{x,t) = 4>{t), this system is equivalent to an isotropic 
harmonic oscillator in two dimensions. This remains the case if we add higher order terms which respect the phase 
symmetry. In supersymmetric models, however, we expect that higher order terms will break the symmetry. In the 
isotropic oscillator analogy, this corresponds to anharmonic terms which break the rotational invariance. With a 
general initial condition, the system will develop some non-zero angular momentum. If the motion is damped, so that 
the amplitude of the oscillations decreases, these rotationally non-invariant terms will become less important with 
time. 

Let us add interactions in the following way, which will closely parallel what happens in the supersymmetric case. 
Include a set of quartic couplings: 

Ci = X\(j)\'^ + e<l>^(l>* +S(p'^ + C.C. (65) 

These interactions clearly violate "B" . For general complex e and S, they also violate CP. In supersymmetric theories, 
as we will shortly see, the couplings X,e,S . . . will be extremely small, 0{Myy/Mp) or 0{M^r/M^^^). 

In order that these tiny couplings lead to an appreciable baryon number, it is necessary that the fields, at some 
stage, were very large. To see how the cosmic evolution of this system can lead to a non-zero baryon number, first note 
that at very early times, when the Hubble constant, _ff 3> m, the mass of the field is irrelevant. It is thus reasonable 
to suppose that at this early time ^ = ^ 0; later we will describe some specific suggestions as to how this might 
come about. How does the field then evolve? First ignore the quartic interactions. In a gravitational background, the 
equation of motion for the field is 

BV 

Dl^+ — = 0, (66) 

where is the covariant derivative. For a spatially homogeneous field, (f>(t), in a Robertson- Walker background, this 
becomes 

4> + 3H^+^=Q. (67) 

At very early times, H ^ m, and so the system is highly overdamped and essentially frozen at (po- At this point, 
B = {). However, once the universe has aged enough that H ^ m, (j) begins to oscillate. Substituting H = or 
iJ = ^ for the radiation and matter dominated eras, respectively, one finds that 

{ , f^°3/2 sin(mt) (radiation) 
f \, , (68) 

^^sm(mt) (matter). 



22 



In either case, the energy behaves, in terms of the scale factor, R{t), as 

E^m'c^li^r (69) 

i.e. it decreases like i?^, as would the energy of pressureless dust. One can think of this oscillating field as a coherent 
state of (/) particles with p^O. 

Now let's consider the effects of the quartic couplings. Since the field amplitude damps with time, their significance 
will decrease with time. Suppose, initially, that = is real. Then the imaginary part of ^ satisfies, in the 
approximation that e and S are small, 

4>^ + 3H^, + m^(|)^^lm{e + S)(f>l (70) 

For large times, the right hand falls as whereas the left hand side falls off only as t~^^'^. As a result, just as in 

our mechanical analogy, baryon number (angular momentum) violation becomes negligible. The equation goes over 
to the free equation, with a solution of the form 

lm{e + d)(j>l . , e \ / . X , Im(f + ())fr' . , e \ / \ 
(pi = ttr — — . „ , . sm{mt + dri (radiation), (pi = am 5 — — sin(mi + Om.) (matter), (71) 

in the radiation and matter dominated cases, respectively. The constants 5m, 5^, am and a^. can easily be obtained 
numerically, and are of order unity: 

ttr = 0.85 am = 0.85 5r = -0.91 5m = 1-54. (72) 

But now we have a non-zero baryon number; substituting in the expression for the current, 

Ub = 2arlm(e + 5) — ° „ sm(5r + 7r/8) (radiation) ub = 2amlm(e + 5) — ° „ siii(5m) (matter). (73) 
mymty mymty 

Two features of these results should bo noted. First, if e and 5 vanish, ub vanishes. If they are real, and (po is 
real, ub vanishes.lt is remarkable that the lagrangian parameters can be real, and yet (po can be complex, still giving 
rise to a net baryon number. We will discuss plausible initial values for the fields later, after we have discussed 
supersymmetry breaking in the early universe. Finally, we should point out that, as expected, ns is conserved at late 
times. 

This mechanism for generating baryon number could be considered without supersymmetry. In that case, it begs 
several questions: 

• What are the scalar fields carrying baryon number? 

• Why are the (p^ terms so small? 

• How are the scalars in the condensate converted to more familiar particles? 

In the context of supersymmetry, there is a natural answer to each of these questions. First, as we have stressed, 
there arc scalar fields carrying baryon and Icpton number. As we will sec, in the limit that supersymmetry is unbroken, 
there are typically directions in the field space in which the quartic terms in the potential vanish. Finally, the scalar 
quarks and leptons will be able to decay (in a baryon and lepton number conserving fashion) to ordinary quarks. 



C. Flat Directions and Baryogenesis 

To discuss the problem of baryon number generation, we first want to examine the theory in a limit in which we 

ignore the soft SUSY-breaking terms. After all, at very early times, H ^ Mw, and these terms arc irrelevant. We 
want to ask whether in a model like the MSSM, some fields can have large vev's, i.e. whether there are directions in 
the field space for which the potential vanishes. Before considering the full MSSM, it is again helpful to consider a 
simpler model, in this case a theory with gauge group U{\), and two chiral fields, </>"'" and (p~ with opposite charge. 
We take the superpotential simply to vanish. In this case the potential is 



v = \d^ d = 9{cp+*(P+ -r*r) 



(74) 
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But D, and the potential, vanish if 0"*" = = v. It is not difRcult to work out the spectrum in a vacuum of non-zero 
V. One finds that there is one massless chiral field, and a massive vector field containing a massive gauge boson, a 
massive Dirac field, and a massive scalar. 

Consider, now, a somewhat more elaborate example. Let us take the MSSM and give expectation values to the 
Higgs and the slepton fields of eqn. (|46(l : 

The F term vanishes in this direction, since the potentially problematic H^L term in the superpotential is absent 
by R parity (the other possible contributions vanish because Q — Hy = 0). It is easy to see that the Z?-term for 
hypercharge vanishes, 

DY = g'^{\H^\^~\m^O. (76) 

To sec that the D terms for SU{2) vanishes, one can work directly with the Pauli matrices, or use, instead, the 
following device which works for a general SU{N) group. Just as one defines a matrix-valued gauge field. 



one defines 



Then, using the SU{N) identity. 



{A.y^^AliT^r^, (77) 

{DYj = D^iry^. (78) 



(r'^)v(T")^-^V'--^'5V^; (79) 

the contribution to [DY j from a field, 0, in the fundamental representation is simply 

(i^)^=0"</>,-l|</.p5i. (80) 

In the present case, this becomes 



2' ' \0 2^ 

What is particularly interesting about this direction is that the field carries a lepton number. As we have seen, 
producing a lepton number is for all intents and purposes like producing a baryon number. 

Non-renormalizable, higher dimension terms, with more fields, can lift the flat direction. For example, the quartic 
term in the superpotential: 

C^ = ^{H^Lf (82) 
respects all of the gauge symmetries and is invariant under i?-parity. It gives rise to a potential 

v.. = 5 (83) 

where $ is the superfield whose vev parameterizes the flat direction. 

There are many more fiat directions, and many of these do carry baryon or lepton number.^ A fiat direction with 
both baryon and lepton number excited is the following; 

First generation : Ql = b U2 = a L2 — b Second : Ji = -I- |ap Third : dz — a. (84) 



The flat directions in the MSSM have been cataloged by iGherghetta. Kolda and Martiri <1996l) . 
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(On Q, the upper index is a color index, the lower index an SU{2) index, and we have suppressed the generation 
indices). 

To check that this is indeed a flat direction, consider first the D terms for the various groups. Using our earlier 
matrix construction, we have: 

/|6p-(|ap + |6|2) \ . 2 n X 

SU{3) : { -|a|2 - trace = S'L/(2) : M ' , ,3 - trace = (85) 



-|o 



2 



|6p 



The F terms also vanish: 



U{1) : l\bf - ^\af + ^\b\' + ^|a|2 + ^\af - \bf = 0. (86) 

-_ = r"V^i' = and —^^'^Q'^dP^Q (87) 
oHu oHd 

The first follows since the u fields have their expectation values in different "color slots" than the Q fields. The second 
is automatically satisfied since the d and Q fields have expectation values in different generations, and these Yukawa 
couplings don't mix generations. 

Higher dimension operators again can lift this fiat direction. In this case the leading term is: 

CT^^lQ^^L^u^S'd']. (88) 

Here the superscripts denote flavor. We have suppressed color and SU{2) indices, but the braces indicate sets of 
fields which are contracted in SU{3) and SU{2) invariant ways. In addition to being completely gauge invariant, this 
operator is invariant under ordinary i?-parity. (There are lower dimension operators, including operators of dimension 
4, which violate i?-parity). It gives rise to a term in the potential: 

$10 

^lift - -J^- (89) 
Here $ refers in a generic way the fields whose vev's parameterize the fiat directions (a, 5). 



D. Evolution of the Condensate 



For the cosmologies we wish to consider, despite the powers of jj, these operators are quite important. During 
inflation, for example, such operators can determine the initial value of the fleld, (here $ denotes in a generic way 
the flclds which parameterize the fiat directions). 



1. Supersymmetry Breaking in the Early Universe 

We have indicated that higher dimension, supersymmetric operators give rise to potentials in the fiat directions. To 
fully understand the behavior of the fields in the early universe, we need to consider supersymmetry breaking, which 
gives rise to additional potential terms. 

We have indicated, in equation (|49|l . the sorts of supersymmetry-breaking terms which we expect in supersymmetric 
theories. In the early universe, we expect supersymmetry is much more badly broken. For example, during inflation, 
the non-zero energy density (cosmological constant) breaks supersymmetry. Suppose that / is the inflaton fleld, and 
that the inflaton potential arises because of a non-zero value of the auxiliary fleld for /, F/ — (see eqn. ()42|l ). F/ is 
an order parameter for supersymmetry breaking as are the auxiliary flelds for any fl eld; this quantity is roughly constan t 
during inflation. So, during inflation, supersymmetry is broken by a large amount ijPine. Randall, and ThomasLlTool . 
Not surprisingly, as a result, there can be an appreciable supersymmetry-breaking potential for These contributions 
to the potential have the form:^ 

Vh = H^'P^fi^^/M^). (90) 



When supersymmetric theories are coupled to gravity, there are corrections to eqn. I41i . It also makes no sense to restrict the 
lagrangian to be renormalizable. Assuming that non-renormalizable couplings scale with Mp, leads to eqn. 1901 . as explained in 
iPine. Randall, and Thomasl IToQ^) . 
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It is perfectly possible for the second derivative of the potential near the origin to be negative. In this case, writing 
our higher dimension term as: 

= (91) 

the potential takes the form 

y=-if2|$|2 + _l_|$|2n+4^ (92) 

The minimum of the potential then lies at: 

" {ji) ■ 

More generally, one can see that the higher the dimension of the operator which raises the flat direction, the larger the 
starting value of the field - and the larger the ultimate value of the baryon number. Typically, there is plenty of time 
for the field to find its minimum during inflation. After inflation, H decreases, and the field $ evolves adiabatically, 
oscillating slowly about the local minimum for some time. 

Our examples illustrate that in models with R-parity, the value of n, and hence the size of the initial field, can vary 
appreciably. With further symmetries, it is possible that n is larger, and even that all operators which might lift the 
flat direction are forbidden l)Dine. Randall, and Thoma s. "1996V For the rest of this section we will continue to assume 
that the flat directions are lifted by terms in the supcrpotcntial. If they are not, the required analysis is different, 
since the lifting of the flat direction is entirely associated with supersymmetry breaking. 



2. Appearance of The Baryon Number 

The term in the potential, does not break either baryon number or CP. In most models, it turns out that 

the leading sources of B and CP violation come from supersymmetry-breaking terms associated with Fj. These have 
the form^ 

am^i/^W + hHW. (94) 

Here a and h are complex^ dimensionless constants. The relative phase in these two terms, 5, violates CP. This is 
crucial; if the two terms carry the same phase, then the phase can be eliminated by a field redefinition, and we have 



to look elsewhere for possible CP-violating effects. Examining equations H82I88|) . one sees that the term proportional 
to W violates B and/or L. In following the evolution of the field <&, the important era occurs when H ~ At 
this point, the phase misalignment of the two terms, along with the B violating coupling, lead to the appearance of 
a baryon number. From the equations of motion, the equation for the time rate of change of the baryon number is 



duB _ sin((5)m3/2 
dt ~ ^ ' 

Assuming that the relevant time is i?^^, one is lead to the estimate (supported by numerical studies) 



(95) 



^B = ^sin(5)$»+3. (96) 



Here, $o is determined by H ^ ^713/2, i.e. <i)2"+2 — m'^^^A'P". 



E. The Fate of the Condensate 

Of course, we don't live in a universe dominated by a coherent scalar field. In this section, we consider the fate of 
a homogeneous condensate, ignoring possible inhomogeneities. The following sections will deal with inhomogeneities, 
and the interesting array of phenomena to which they might give rise. 



Again, these arise from non-rcnormalizable terms in the effective action, <Dine. Randall and Thomad IToQ^) . 
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We have seen that a coherent field can be thought of as a collection of zero momentum particles. These particles 
are long-lived, since the particles to which they couple gain large mass in the flat direction. Were there no ambient 
plasma or other fields, the condensate would eventually decay. However, there are a number of effects which lead 
the condensate to disappear more rapidly, or to produce stable remnants. Precisely what arc the most important 
mechanisms depend on a number of factors. Perhaps most important is the rate of expansion, and the dominant form 
of energy during this epoch. The amplitude of oscillations is also important. 

It is impossible to survey all possibilities; indeed, it is likely that all of the possibilities have not yet been imagined. 
Instead, we will adopt the picture for infiation described in the previous section. The features of this picture are true 
of many models of inflation, but by no means all. We will suppose that the energy scale of inflation is i? ~ 10^^ GeV. 
We assume that inflation is due to a field, the inflaton I. The amplitude of the infiaton, just after inflation, we will 
take to be of order M ^ 10^^ GeV (the so-called reduced Planck mass). Correspondingly, we will take the mass of 
the inflaton to be mj = 10"'^^ GeV (so that rnjMp « E'^). Correspondingly, the Hubble constant during inflation is of 
order Hi « E^/Mp « lO^^ Q^y. 

After inflation ends, the inflaton oscillates about the minimum of its potential, much like the field $, until it decays. 
We will suppose that the inflaton couples to ordinary particles with a rate suppressed by a single power of the Planck 
mass. Dimensional analysis then gives for the rough value of the inflaton lifetime: 

3 

T/ = ^ ^ 1 GcV. (97) 

The reheating temperature can the n be obtained by equatin g the energy density at time {H ^ T; p — 3H^AP to the 
energy density of the final plasma: ijKolb and Turnerl 11990(1 

Tr ^ T{t ^ rj^) - (r/Mp)i/2 - lO^GeV (98) 

The decay of the inflaton is actual ly not sudden, but leads to a gradual reheating of the imiverse, as described, for 
example, in ijKolb and Turneil Il990|) . As a function of time (H): 

T « {TlH{t)Mpf'^. (99) 

As for the field our basic assumption is that during inflation, it obtains a large value, in accord with equation 
(|93|l . When inflation ends, the inflaton, by assumption still dominates the energy density for a time, oscillating about 
its minimum; the universe is matter dominated during this period. The fleld now oscillates about a time-dependent 
minimum, given by equation (|93|l . The minimum decreases in value with time, dropping to zero when H ^ 
During this evolution, a baryon number develops classically. This number is frozen once H ~ 

Eventually the condensate will decay, through a variety of processes. As we have stressed, the condensate can be 
thought of as a coherent state of <& particles. These particles - linear combinations of the squark and slepton fields 
- are unstable and will decay. However, for H < the lifetimes of these particles are much longer than in the 

absence of the condensate. The reason is that the fields to which <& couples have mass of order $, and $ is large. In 
most cases, the most important process which destroys the condensate is what we might call evaporation: particles 
in the ambient thermal bath can scatter off of the particles in the condensate, leaving final states with only ordinary 
particles. 

We can make a crude estimate for the reaction rate as follows. Because the particles which couple directly to <& are 
heavy, interactions of $ with light particles must involve loops. So we include a loop factor in the amplitude, of order 
the weak coupling squared. Because of the large masses, the amplitude is suppressed by $. Finally, we need to 
square and multiply by the thermal density of scattered particles. This gives: 

~ ^l^^iTRHMf/*. (100) 

The condensate will evaporate when this quantity is of order H. Since we know the time dependence of $, this allows 
us to solve for this time. One finds that equality occurs, in the case n — 1, for Hj ~ 10^ — 10'^ GeV. For n > 1, it 
occurs significantly later (for n < 4, it occurs before the decay of the inflaton; for n > 4, a slightly different analysis 
is required than that which follows). In other words, for the case n ~ 1, the condensate evaporates shortly after the 
baryon number is created (but for more complications, see below), but for larger n, it evaporates significantly later. 

The expansion of the universe is unaffected by the condensate as long as the energy density in the condensate, 
p$ ~ m|,$^, is much smaller than that of the inflaton, pi ~ H^M^. Assuming that mg, ^ 7713/2 ~ 0.1 — 1 TeV, a 
typical supersymmetry breaking scale, one can estimate the ratio of the two densities at the time when H ^ SiS 



I -m , \ 2/(n+l) 



PI V M 



(101) 
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We are now in a position to calculate the baryon to photon ratio in this model. Given our estimate of the inflaton 
lifetime, the coherent motion of the inflaton still dominates the energy density when the condensate evaporates. The 
baryon number is just the <I> density just before evaporation divided by the $ mass (assumed of order 7713/2), while 
the inflaton number is pi /Mj. So the baryon to inflaton ratio follows from eqn. I|1U1|) . With the assumption that the 
inflaton energy density is converted to radiation at the reheating temperature, T/?, we obtain: 

77b ^ ns ^ nB_TR_p^ ^ 10"^" ( \ ( 

(pi/Tr) n^mq,pi VlO^GeV/ V"^3/2 

Clearly the precise result depends on factors beyond those indicated here explicitly, such as the precise mass of the 
$ particle(s). But as a rough estimate, it is rather robust. For 71 = 1, it is in precisely the right range to explain 
the observed baryon asymmetry. For larger 71, it can be significantly larger. While this may seem disturbing, it is 
potentially a significant virtue. Many supersymmetric models lead to creation of entropy at late times. For example, 
in string theory one expe cts the existence of oth er "light" {m « 7713/2) fields, known as "moduli." These fields lead to 
cosmological difficulties llCouehlan et a/.l . Il983|) . unless, when they decay, they reheat the universe to temperatures 
of order 10 MeV, after which nucleosynthesis can occur. These decays produce a huge amount of entropy, typically 
increasing the energy of the universe by a factor of 10^. The baryon density is diluted by a corresponding factor. So 
in such a picture, it is necessary that the baryon number, prior to the moduli decay, should be of order 10^'^. This is 
not the only cosmological model which requires such a large baryon number density. 

There are many issues in the evolution of the condensate which we have not touched upon. One of the most ser ious 
is related to interactions with the thermal bath ( AUah verdi. Campbell, and Ellij . I2OOOI : lAnisimov and Dind.[2000() . In 
the case 71 = 1, $0 is not so large, and, while the particles which $ couples to get mass of order 4>o, they may be in 
thermal equilibrium. In this case, the <& particles decay much earlier. This typically leads to significant suppression 
of the asymmetry, and the viability of the AD mechanism depends on the precise values of the parameters. 

Overall, then, there is a broad range of parameters for which the AD mechanism can generate a value for ^ equal 
to or larger than that observed. This baryon number is generated long after inflation, so inflationary reheating does 
not provide any significant constraint. It can be large, allowing for processes which might generate entropy rather 
late. 



F. Inhomogeneities and the Condensate 

We have so far assumed that the condensate is homogeneous. But, as we will now show, under some circumstances 
the condensate is unstable to fragmentation. This appears to be related to another feature of theories with scalars: 
the possible existence of non-topological solitons. These can alter the picture of baryon number generation, and could 
conceivably be a dark matter candidate. This is the subject of this section. 

1. Stability and fragmentation 

To analyze the stability of the condensate l)Kusenko and Shaposhnikovlll99^ . we write the complex field (p — pe*^^ 
in terms of its radial component and a phase, both real functions of space-time. We are interested in the evolution 
of the scalar field in the small- VEV domain, where the baryon number violating processes are suppressed, and we 
will assume that the scalar potential preserves the U(l) symmetry: U{lp) = U{p), where U{p) may depend on time 
explicitly. The classical equations of motion in the spherically symmetric metric c?s^ dt^ — a{t)'^dr'^ are 

n + iHh- ^—An+'^ti- ^^{d^^){d'p) = 0, (103) 

a-^it) p O' \t)p 

p + 3Hp- - n^p + -L_{d,nfp + {du/dp) - 0, (104) 

where dots denote time derivatives, and the space coordinates are labeled by the Latin indices that run from 1 to 3. 
The Hubble constant, again, is H = a/a, where a{t) is the scale factor; it is equal to t~^^^ or t~^/^ for a matter or 
radiation dominated universe, respectively. 

Prom the equations of motion (|103|l and (|1U4|I . one can derive the equations for small perturbations Sil and Sp: 



(102) 
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FIG. 5 The charge density per comoving volume in (1+1) dimensions for a sample potential analyzed numerically during the 
period when the spatially homogeneous condensate breaks up into high- and low-density domains. Two domains with high 
charge density are expected to form Q-balls. 



(517 + 3i7(^0) - ^r--A((5r2) + — (517) + — {8p)~Ar^P = 0, (105) 
a^(t) p p p'^ 

5p + 3H(Sp)- ^^A(5p)-2pn(Sh) + U"Sp-n^6p = 0. (106) 

To examine the stability of a homogeneous solution ip{x,t) — ip{t) = p(i)e*^(*', let us consider a perturbation 
Sp,Sn oc QS(t)-ikx \ooY for growing modes, Re a > 0, where a = dS/dt. The value of k is the spectral index 
in the comoving frame and is red-shifted with respect to the physical wavenumber k — k/a(t) in the expanding 
background. Of course, if an instability develops, the linear approximation is no longer valid. However, we assume 
that the wavelength of the fastest-growing mode sets the scale for the high and low density domains that eventually 
evolve into Q-balls. This assumption can be verified post factum by comparison with a numerical solution of the 
corresponding partial differential equations (|103|l and (|104|l . where both large and small perturbations are taken into 
account. 

The dispersion relation follows from the equations of motion: 



a + Ana H — H a 





P 

a 




P. 



a = 0. 



(107) 



If (f2^ — U"{p)) > 0, there is a band of growing modes that lies between the two zeros of a(fc), < fc < fcmax, where 



km.At) = a{t)^h^^U"{p). (108) 

This simple linear analysis shows that when the condensate is "overloaded" with charge, that is when u){t) — Cl 
is larger than the second derivative of the potential, an instability sets in. Depending on how A:max(i), defined by 
relation H108|l . varies with time, the modes in the bands of instabil i ty may or may not have time to develop fully. 

Numerical analyses iEn avist el oi.L 120011: iKasuva and Kawasakii l2000albl l200ll) , which can trace the evolution of 
unstable modes beyond the linear regime, have shown that fragmentation of the condensate is a generic phenomenon. 
Numerically one can also study the stability of rapidly changing solutions, hence relaxing the adiabaticity condition 
assumed above. This aspect is relevant to the cases where the baryon number density is small and the radial 
component of the condensate, p(i), exhibits an oscillatory behavior changing significantly on small time scales. An 
interesting feat ure of this non-adiab atic regime is that both baryon and anti-baryon lumps may form as a result of 
fragmentation l)Enavist et a/.U200l[l . 
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2. Lumps of scalar condensate: Q-Balls 



Perhaps the most famihar sohton solutions of non-linear field theories, such as magnetic monopoles and vortices, 
can be uncovered by to pological ar guments. However, fi eld theories with scalar fields often admit non-topological 
solitons lIColeinanL Il98,5: Friedberg. Lee, and SirlinL Il97fit iRosen. 1968.'). Q-balls, which may be stable or may decay 
into fermions ijCohen. Coleman. Georgi. and Manohailll986(l . Q-balls appear when a complex scalar field ip carries 
a conserved charge with respect to some global U(l) symmetry. In supersymmetric generalizations of the Standard 
Model, squarks and sleptons, which carry the conserved baryon and lepton numbers, can form Q-balls. 

Let us consider a field theory with a scalar potential Ui^p) which has a global minimum ?7(0) = at 1^9 = 0. Let 
V{}p) have an unbroken global U(l) symmetry at the global minimum: p — > exp{i0}i^. We will look for solutions of 
the classical equations by minimizing the energy 



(109) 



subject to the constraint that the configuration has a definite charge, Q, 

Q ^ -^f If* dt pd^x 



(110) 



To describe the essential features of Q-balls in a simple way, we will, following IColemanI l)l985j) . use a thin wall 
ansatz for the Q-ball 



(/.(x,t) = e^"'0(a;), (111) 

where 

' 0, > R 



= \ : (112) 

(for the real solution, the field varies rapidly between the two regions, changing on a scale of order the Compton 
wavelength of the 4> particle). 

Assuming that Q is large, let us neglect the gradient terms (relevant only for the wall energy). Then the global 
charge and the energy of the field configuration I|111I112|I are given by 



uJ^oV. (113) 



where V = {4/3)ttR^, and 



E=^Lo'cj,lV + U{MV^ly^^+VU{M (114) 
We now minimize E with respect to V, obtaining 



It remains to minimize the energy with respect to variations oi (fio- A non-trivial minimum exists as long as 



U{ip) /ip'^ ~ min, for (p = ipa > 0; (116) 

if this condition is satisfied, a Q-ball solution exists. 

So far we have assumed a particular ansatz, neglected the gradient terms, etc. These assumptions can be avoided in 
a slightly more involved derivation (,KusenkOi. .1997a^) using the method of Lagrange multipliers. We want to minimize 
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Q 



1 

2i 



(117) 



where w is a Lagrange multiplier. (It is no accident that we use the same letter, w. The value of the Lagrange 
multiplier at the minimum will turn out to be equal to the time derivative of the phase.) Variations of ip{x,t) and 
those of uj can now be treated independently, the usual advantage of the Lagrange method. 
One can re-write equation H117|l as 
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/a 



where 



-|V¥.p + ;7^(^) 



(118) 



2 2 
■UJ (p . 



(119) 



We are looking for a solution that extremizes while all the physical quantities, including the energy, E, are 
time-independent. Only the first term in equation H118|) appears to depend on time explicitly, but it is positive definite 
and, hence, it should vanish at the minimum. To minimize this contribution to the energy, one must choose, therefore. 



(^(x,t) 



V(2^), 



(120) 



where (p{x) is real and independent of time. We have thus derived equation (|111|) . For this solution, equation H110|) 
yields 



u) / (p'^{x) d^x 



(121) 



It remains to find an extremum of the functional 



-\Vip{x)\' + UUv'ix)) 



ujQ, 



(122) 



with respect to lu and the variations of (p{x) independently. We can first minimize £^ for a fixed u, while varying the 
shape of (p{x). If this were an actual potential for a scalar field in three dimensions, one would have the possibility 
of tunneling between the zero energy configuration at the origin and possible lower energy configurations at non-zero 
(j) (Fig. 6). Tunneling, in the semiclassical approximation, is described by the bounce , (pi^{x), the solution of the 
classical equations which asymptotes to the "false vacuum" at the origin The first term in equation 

(|122|l is then nothing but the three-dimensional Euclidean action 5*3 [(^^^(x)] of this bounce solution. This is a very 
useful correspondence. In particular, the condition for the existence of solution is simply a corollary: as long as Uuj{(p) 
has a minimum below zero, the bounce exists, and so does the Q-ball, ip{x,t) = exp{iu!t}(p{x). 

The bounce, and hence the Q-ball, exist if there exists a value of w, for which the potential Ui^{ip) has both a 
local minimum at (p = and a global minimum at some other value of if. This condition can be re-phrased without 
reference to w: a Q-ball solution exists if 



U{ip) I 'P^ = min, for (y9 = > 



(123) 



The corresponding effective potential U^^^^iip)^ where = \/ "^U {ip^) / ip^ , has two d egenerate minima, at = and 
ip = ipo- The existence of the bounce solution (puj{x) for coq < lu < U"{0) follows ijColemanL Il977t IColeman et all 
Il978l) from the fact that tJuif) has a negative global minimum in addition to the local minimum at the origin. 
IColeman et al\ (^23) also showed that the solution is spherically symmetric: pix) = ipir), r ~ . 

The soliton we want to construct is precisely this bounce for the right choice of w, namely that which minimizes 
The last step is to find an extremum of 
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~- large co, small Q 

FIG. 6 Finding a Q-ball is equivalent to finding a bounce that describes tunneling in the potential Uij{ip) — U{ip) — {l/2)Lu'^ip^ . 
The thin- wall approximation is good for large Q (upper dashed line), but breaks down when Q is small and, therefore, u is 
almost as large as the mass term at the origin. In the latter case (lower dashed line), the "escape point", <fi{0) is close to the 
zero of the potential and is far from the global minimum. 



£u, = Ss[^^{x)] + ojQ (124) 

with respect to uj. One can prove the existence of such an extremum jKusenkol Il997a|) . Finally, the soliton is of the 
form (|120|l . with uj that minimizes £ui in eq. H124II . 

Having obtained the solution, one can compute its energy (mass). For a finite ipo in eq. ( 1123(1 . in the limit of large 
Q, the Q-ball has a thin wall, and its mass is given by 



M{Q) = LJnQ (125) 

Supersymmetric generalizations of the Standard Model have scalar potentials with flat directions lifted only by 
supersymmetry breaking terms. Q-balls may form with a light scalar field ip that corresponds to that flat direction. 
If the potential V{(p) = (x^ = const for large then the minimum in eq. I|123|) is achieved for (/Jq = oo. In this case, 

M{Q) = ^xQ'"\ (126) 

More generally, if the potential grows slower than (/s^, i.e. V{(.p) oc 93^, p < 2, condition ((123|1 is not satisfied at any 
finite value of 0Oi and 



M{Q) - ^Q(3-p/2)/(4-p)^ (^27) 

It is important in what follows that the mass per unit charge is not a constant, but is a decreasing function of the 
total global charge Q. There is a simple reason why the soliton mass is not proportional to Q. Since U(ip)/f'^ has no 
minimum, the scalar VEV can extend as far as the derivative terms allow it. When the next unit of charge is added, 
the Q-ball increases in size, which allows the scalar VEV to increase as well. Hence, the larger the charge, the greater 
is the VEV, and the smaller is energy per unit charge. 



3. AD Q-balls 

The condition H123|) suggests that slowly growing potentials, of the sort that arise in the flat directions of the MSSM, 
are a likely place to find Q-balls. Q-balls can develop along flat directions that carry non-zero baryon number, lepton 
number, or both. Each flat direction can be parameterized by a gauge-invariant field, carrying these global quantum 
numbers. So the discussion of gauge singlet fields of the previous section also applies to baryonic and leptonic Q-balls 
in the MSSM. This statement may seem surprising, since all scalar baryons in the MSSM transform non-trivially 
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under the gauge group. Although scalars with gauge interactions c an also make Q-balls (|L ee et al l Il989l), in th e 
case of the MSSM the color structure of large Q-balls is rather simple (jKusenko. Shaposhnikov. and Tinvakovlll998|) . 
If a Q-ball VEV points along a flat direction, its scalar constituents form a colorless combination (otherwise, that 
direction would not be flat because of non- vanishing D-terms). 

In one proposal for the origin of supersymmetry breaking, "gauge-mediated" supersymmetry breaking, breaking, 
the flat directions are lifted by potentials which grow quadratically for small values of the fields, and then level off to 
a logarithmic plateau at larger 0. Q-balls in such a potential have masses given by eq. H12t)|) . In another proposal, 
"gravity-mediated" scenarios, the potentials which arise from supersymmetry breaking grow roughly quadratically 
even for very large VEV. Whether Q-balls exist is thus a detailed, model-dependent question. Q-balls in these 
potentials have masses proportional to the first power of Q. 

By con struction, Q-balls are stable with respect to decay into scalars. However, they can decay by emitting 
fermions ijCohen. Coleman . Georgi. and Manohar. 1986). If the Q-ba ll has zero baryon number, it can decay by 
emitting light neutrinos jCoheiTColemair^eorgi. and Manohail Il986l) . 

However, if a baryonic Q-ball ("B-ball") develops along a flat direction, it can also be stable with respect to decay 
into fermions. Stability requires that the baryon number be large enough. A Q-ball with baryon number Qb and 
mass M{Qb) is stable if its mass is below the mass oi Qb separated baryons. For a Q-ball in a flat potential of height 
Ms^ the mass per unit baryon number 

~ MsQ-'" (128) 

Models of gauge- mediated supersymmetry breaking produce flat potentials with A/5 1 — 10 TeV. If the mass per 
baryon number is less than the proton mass, Wp, then the Q-ball is entirely stable because it does not have enough 
energy to decay into a collection of nucleons with the same baryon number. This condition translates into a lower 
bound on Qb'- 



4. Dark matter in the form of stable B-balls 

Stable Q-balls that form from an AD condensate are a viable candidate for dark matter. Even if they are unstable, 
their decay can produce neutralinos at late times, when these neutralinos are out of equilibrium. One way or another, 
some dark matter can arise from AD baryogenesis. 

Moreover, since both the ordinary matter and the dark mat ter have the same origin i n the AD scenario, one can try 
to explai n why their am ounts in the univ erse are fair ly close (IBaneriee and JedamzlH l200ft lEnavist and McDonaldL 
[1999; Fui ii and Yanagid a. 2002at iLainea nd Shaposhi iikovl Il998|) . 

Since the MSSM with gauge-mediated supersymmetry breaking contains stable objects, baryonic Q-balls, it is 
natural to ask whether they can constitute the dark matter. Stable Q-balls can be copiously produced in the course of 
AD baryogenesis. Of the dark-matter candidates that have been considered, most were weakly interacting particles, 
and for good reason. If dark-matter particles have strong [relative to their mass) interactions with matter, these 
interactions might facilitate their loss of momentum and angular momentum, forcing them into the galactic disks, 
along with ordinary matter. But astronomical observations show that the dark matter forms spherical halos about 
galaxies, not disks. 

Made of squarks, these Q-balls can interact strongly with ordinary matter, via QCD. However, if they are as heavy as 
the calculations show they are, the strong interactions are not strong enough to force dark-matter Q-balls to settle into 
the galactic disks. Analyses of Q-ball formation and partial evaporation allow one to relate the amounts of ordinary 
matter and dark matter. The observed ratio corresponds to Q-balls w ith baryon number of about 10^^^^, which is in 
agreement with the expected Q-balls size from numeri cal simulation s ( Kasu va and Kawasakil l2000al l200l|) , as well as 
with the current experimental bounds summarized bv lArafune et aL (:2000,) . see Fig.EI A B-ball with baryon number 
10^^ is so heavy that it could pass through ordinary stars with only a small change in its velocity {{5v/v) ~ 10~^). 
Hence, despite the strong interactions, B-balls make a good dark-matter candidate. 

The unusually large mass for a dark matter candidate means that the fluxes are very small. The experimental 
detection of this form of dark matter requires a large detector size (c/. Fig. [7||. 

Since B-balls have lower mass-to-baryon ratio than ordinary nuclear matter, interactions of B-balls with ordi- 
nary matter result in numerous "proton decay" events (Kuscnko et ai, 1998a). Hence a Q-ball passing through 
a detector would produce a spectacular signature. The flux, however, is very low. Hence, the strongest limits 
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come from the largest detec tors, e.g. Super-Kamiokande (see Fig. EJ. Some astrophysical bounds have been consid- 
ered (Kuscnko et al, 1998c), but they do not yield very strong constraints. In addition to the existing limits discussed 
in (jArafune et al... .2000.') . future experiments, such as ANTARES, Ice Cube, etc, may be able to detect dark matter 
B-balls, or rule out the values of Q that correspond to the correct amount of dark matter. 

We note that, although Q-balls are always present in the spectrum of any SUSY extension of the Standard Model, 
their production in the early universe requires the formation of an Affleck- Dine condensate followed bv its fragme nta- 
tion. Stable Q-balls are too large to form in thermal plasma by accretion l|Griest and Kolbl Il989l: lPostmall2002(l . In 
this sense, an observation of stable dark-matter Q-balls would be evidence of the Affleck-Dine process having taken 
place. 



5. Dark matter from unstable B-balls 

If supersymmetry breaking is mediated by gravity, Q-balls are not stable as they can decay into fcrmions. However, 
Q-ball decay into fcrmions is a slow process because the fermions quickly fill up the Fermi sea inside the Q-ball, 
and further decays are limited by the rate of fermion evaporation through the surface. The rate of Q-ball decay is, 
therefore, suppressed by the surface-to- volume ratio ( Cohen, Coleman. Georgi. and Manphar,, 198Q) as compared to 
that of free scalar particles. In a typical model, unstable baryonic Q-balls from the Affleck-Dine condensate decay 
when the temperature is as low as a GeV. 

The lightest supersymmetric particles (LSP) are among the decay products of Q-balls. B-balls can decay and 
produce dark matter in the form of neut ralinos at a time when they are out of equilibrium ijEnovist and McDonaldL 
Il998l IT999t iFuiii and HamaguchiL l2002al) . This presents another possibility for producing dark matter from the AD 
condensate and relating its abundance to that of ordinary matter. The requirement tha t neutralinos not overdose 
the universe constrains the parameter space of the MSSM ijFuiii and Hamag uchi. '2002a'' b|). 

If the LSPs don't annihilate, the ratio of ordinary matter to dark matter is simply lEnavist and McDonaldL Il999|) 

a„attcr/f^LSP ^ /"^ ( ( — ) (130) 



J \ n 



where / is the fraction of the cond ensate trapped in Q-balls If / ^ 10^^, this ratio is acceptable. 

However, numerical simulations ("Kasuva and Kawasakil l2000afl and some analytical calculations l|McDonaldL l200ll) 
indicate that, in a wide class of AD models, practically all the baryon number may be trapped in Q-balls, that 
is / ~ 1. If that's the c ase, the LSP would overdose the universe, according to eq. H13U|) . A solution, proposed 
by Fuiii and HamaeuchJ l)2002allj ). is to eliminate the unwanted overdensity of neutralinos by using an LSP with a 



higher annihilation cross section. The LSP in the MSSM is an admixture of several neutral fermions. Depending 
on the parameters in the mass matrix, determined largely by the soft SUSY breaking terms, the LSP can be closely 
aligned with Bino (one of the SUSY partners of the SU(2)xU(l) gauge bosons), Higgsino (the fermion counterpart of 
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FIG. 8 The allowed range of parameters for non-thermal LSP dark matter is very different from that in the standard freeze-out 
case. The lightly shaded region abov e the solid an d the dashed linos is allowed for non-thermal LSP dark matter in the minimal 
supergravity model with tan /3 = 40 l|Fuiii and HamaeuchL ■2002a .b,'l . 
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FIG. 9 The fate of the AD condensate 



the Higgs boson), or with one of the other weak eigenstates. The trad itional, freeze-out scenari o for LSP production 
favors th e Bino-Kke LSP ijJungman. Kamionkowski. and Griestlll996|) . However, according to iFuiii and Hamaeuchil 
l|2002alb() . SUSY dark matter produced from the Af5eck-D ine process has to be in the form of a Higgsino-Hke LSP. 
In this case, the ratio of matter densities is ((Fuiii and Yanagida. 2002a'l 

anattcr/f^LSP - 10^-' ( 7^ | ( — ) ^CP, (131) 

where ^cp ~ 0.1 is the effective CP violating phase of the AD condensate. For a Higgsino-like LSP, which has 
(av) ~ lO^'^^^^^GeV"^, this yields an acceptable result. 

This has important implications for both direct dark matter searches and the collider searches for SUSY. First, 
the parameter space of the MSSM consistent with LSP dark matter is very different, depending on the cosmological 
scenario at work, that is whether t he LSPs froze out of equili brium tfArnowit t and Dutta , 2002) or were produced 
from the evaporation of AD B-balls ijFuiii and Yanagidall2002af) . Second, higgsino and bino LSP's interact differently 
with matter, so the sensitivity of direct dark-matter searches also depends on the type of the LSP. 

If supersymmetry is discovered, one will be able to determine the properties of the LSP experimentally. This will, 
in turn, provide some information on the how the dark-matter SUSY particles could be produced. The discovery of 
a Higgsino-like LSP would be a evidence in favor of Affleck-Dine baryogenesis. This is yet another way in which we 
might be able to establish the true origin of matter-antimatter asymmetry. 
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IV. CONCLUSIONS 

The origin of the matter-antimatter asymmetry is one of the great questions in cosmology. Yet we can obtain only 
limited information about the events which gave rise to the baryon asymmetry by looking at the sky. Filling out 
the picture requires a deeper understanding of fundamental physical law. One elegant possibility, that the Minimal 
Standard Model produced the baryon number near the electroweak scale, is ruled out decisively by the LEP bounds 
on the Higgs mass. This is a bittersweet conclusion: while one has to give up an elegant scenario, this is perhaps the 
strongest evidence yet for physics beyond the Standard Model - a precursor of future discoveries. 

Supersymmetry is widely regarded as a prime candidate for such new physics. Theoretical arguments in favor of su- 
persymmetry are based on the naturalness of the scale hierarchy, the success of coupling unification in supersymmetric 
theories, and the nearly ubiquitous role of supersymmetry in string theory. The upcoming LHC experiments will put 
this hypothesis to a definitive test. If low energy supersymmetry exists, there are several ways in which it might 
play the crucial role in baryogenesis. It could conceivably revive the electroweak baryogenesis scenario. However, the 
phase transition in the MSSM is only slightly stronger than that in the Standard Model; a noticeable improvement 
forces one into a narrow corner of the MSSM parameter space, which may soon be ruled out. 

But supersymmetry opens a completely new and natural avenue for baryogenesis. If inflation took place in the 
early universe, for which we have an increasing body of evidence, then formation of an Affleck-Dine condensate and 
subsequent generation of some baryon asymmetry is natural. In a wide class of models this process produces the 
observed baryon asymmetry. Perhaps more striking is that the process can lead to very large baryon asymmetries. 
This may be important in many cosmological proposals where one produces substantial entropy at late times. 

Finally, the same process can produce dark matter, either in the form of stable SUSY Q-balls, or in the form of 
a thermally or non-thermally produced LSP. There are hints that the relative closeness of matter and dark matter 
densities may find its explanation in the same process as well. If supersymmetry is discovered, given the success of 
inflation theory, the Affleck-Dine scenario will appear quite plausible. 

Other independent indications that the Affleck-Dine process took place in the early universe may come from 
detection of dark matter. One of the great attractions of supersymmetry is that it can naturally account for the 
dark matter. The lightest supersymmetric particle (LSP) is typically stable, and is produced with an abundance 
in a suitable range if the supersymmetry breaking scale is of order lOO's of GeV. Its precise contribution to the 
energy density of the universe depends on its annihilation cross section and mass. A combination of accelerator limits 
and cosmology presently allows for LSP dark matter in a range of parameters. In this range, the LSP, which is an 
admixture of several states, must be principally "Bino" (the supersymmetric partner of the U{1) gauge boson) if it is 
produced in the standard freeze-out scenario. However, if future detection will indicate that the LSP is Higgsino-like 
(i.e. primarily the partner of the Higgs boson), this kind of dark matter could only arise f rom non-thermal productio n 
of the LSP from a fragmented Affleck-Dine condensate l)Fuiii and Hamaeuchl l2002albl : iFuiii and YanaeidaL l2002a)) . 
Therefore, although a standard Bino-like LSP is not inconsistent with the Affleck-Dine scenario, a Higgsino-like LSP 
would provide a strong evidence in its favor. Likewise, a detection of stable baryonic Q-balls would be a definitive 
confirmation that an Affleck-Dine condensate formed in the early universe and fragmented into B-balls. Since stable 
SUSY Q-balls must be large, we know of no other cosmological scenario that could lead to their formation. 

Among other possibilities for baryogenesis, leptogenesis is also quite plausible. The discovery of neutrino mass, 
perhaps associated with a rather low scale of new physics, certainly gives strong support to this possibility. The 
questions of what scales for this physics might be compatible with inflation, and what implications this might have 
for the underlying origin of neutrino mass are extremely important. Some pieces of the picture will be accessible to 
experiment, but many of the relevant parameters, including the relevant CP violation, reside at a very high scale. 
Perhaps, in a compelling theory of neutrino flavor, some of these questions can be pinned down. 

Future experimental searches for supersymmetry, combined with the improving cosmological data on CMBR and 
dark matter, will undoubtedly shed further light on the origin of baryon asymmetry and will provide insight both 
particle physics and cosmology. The study of the baryon asymmetry has already provided a compelling argument for 
new physics, and holds great promise of new and exciting discoveries in the future. 
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